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1.0 ABSTRACT 


In recent years, data have been obtained which show 
that changes in the chemical structure of the binder of a 
propellant can be correlated with changes in the mechanical 
properties of the propellant. It has also been demonstrated 
that these changes to the chemical structure of the binder 
can be detected using infrared spectroscopy. ‘Therefore, the 
changes in the infrared absorotion characteristics of the 
binder are directiy related to changes in the chemical struc- 
ture and, hence, to changes in mechanical properties of the 
propellant. The computer programs described in this manual 
were developed to demonstrate the feasibility of using Fourier 
transform infrared spectroscopy as a non-destructive tool to 
relate changes in binder structure to changes in mechanical 
properties of the propellant and eventually to use thet 
correlation to predict the remaining service life of the 


propellant. 
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2.0 OVERVIEW OF THE FTIS COMPUTER PROGRAMS 


The Fourier Transform {nfrared Spectroscopy (FTIS) 
computer codes described in this manual,vwere written in 
FORTRAN IV - H EXTENDED for use on either CDC 6000 series or 
IBM 360/370 series digital computers. The primary purposes 
of these codes are to reduce infrared spectral data stored 
on a Digilab FTS-10 tape and correlate this spectral informa- 


tion with changes in propellant physical properties. 


In order to accomplish this goal, the necessary two 


codes are: 


I. E490 - Performs the Infrared data reduction and 
statistical anaiysis. 


II. £410 - Generates a master file for Mechanical 
Property data. 
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2.1 ORGANIZATION OF THE E410 Ano E490 COMPUTER COLES 


The illustration below represents a funct‘cnal flowchart 


for both the E410 and E490 programs. 


Logical Sequence in which 
Subroutines are called bv 


E490MAIN. 

H490MALN 
AL ILF* (CbC Version Only). 
NREAD* (CDC Version Only). 


PLWAVE {If input variable 
IPLOT = 1). 


APCALC 


PLWAVE (If input variable 
NMPLOT = 1). 


BASELN 


COBASE (If input ‘ariable 
MULTPK # 0). 


> 
eos 
PKTABL 
zs REGRES (if input variable 
{ END MPDISK = 1). 


£410: 


E410MAIN 


RENAME (If updating or 
accessing an existina 
END ) properties file). 
y 


*Please see description of 
E490MAIN (3.0.1). 
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3.0 £490 - COMPUTER PROGRAM DISCRIPTION 


The following are general descriptions af both the 


main program plus the six (6) svbroutines that comprise 


the E490 computer code. 


3.0.1 MAIN PROGRAM - E490MAIN 


The MAIN program controls the flow of calculations 
for this FTIS spectral arulysis utilizing both the NAME- 
LIST input supplied by the user and a direct linkage 
to the six subroutines. Other key operations performed 


by E490MAIN are described below: 


a). Printing of the input variables supplied 
by the NAMELIST called /FILES/. 


b). Reads the Digilab FTS-i0 data tape. 
An FTS tape file consists of two physical 
records...a header remarks reccerd and a 
spectral data record. From the header 
record, the following vzriables are read: 


IPFILE() - The file number. 

MWORDS () - File size in data words. 

NSPR - Single / dovble precision designato:. 

NEXP - Exponent of two; (used to compute FTLS 
amplitude values). 

Iti - Wave number (frequency) of the first 
point. (Integer portion;. 

II2 - Numeratcr of the decimal fraction for IIl. 

I13 - Wave number of the last point. 
(Integer portion). 

Ir4 - Numerator of the decirai fraction for [I3. 


The spectral data recerd contains the raw spectrai 
data in binary integer form. 
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c). 


d). 


e). 


f). 
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NOTE: Because of the internal differences 
between IBM and CDC ccnputer systems, the 
current method for reading this FTS-10 tape 
on the Air Force /RPL CDC 6600 computer is 
different. On tne CDC system, E490MAIN calls 
two additional subroutines: XFILE and XREAD. 


XFILE : Called at the beginning of E490MAIN, 
XFILE reads the entire FTS-10 input 
tape and copies it to a temporary file. 


XREAD : By accessing the temporary file gen- 
erated by XFILE, this routine will 
read a record from an FTS file. 


Calculation of the wave number (WAVENO (i) ) 
and amplitude (AMPLTD(i)) for each spectral 
data point i within a file. 


Tabulation of these wave numbers and amplitudes. 


Normalization of these amplitvdes, utilizing 
information computed in subroutine APCALC. 


Tabulation of thece normalized amplitudes and 
their corresponding wave numbers. 
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3.0.2 SUBROUTINE APCALC 


This subroutine calculates variables necessary for the 
normalization of the amplitude values (AMPLTD(i)) in all 
spectral files. The normalization of the array AMPLTD(i) 
is computed relative to the peak with the maximum amplitude 
nearest WAVNOR, a wave number value input by the user. For 
each individual FTIS file, the following normalization para- 


meters are computed in APCALC then returned to the E490MAIN 


routine: 
Variable Name Definition 
a) APMAX Maximum amplitude nearest WAVNOR. 
b) WE Wave number at APMAX, 
c) D2 Calculated amplitude at WE along the 
baseline of the peak. 
ad) AP2 Height of the peak located at WE. 
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3.0.3 SUBROUTINE BASELN 


Subroutine BASELN searches for and records up to a 


maximum of fifty ‘/50) peak heights within each individual 


spectrum. Due to the physical model of this project, the 
search for peaks is done only over certain wave number 
intervals, which are defined below: 

I. 3200. = Wave Numbers = 2700. 


II. 1800. = Wave Numbers = 700. 


NOTE: For files with a last Wave Number greater than 700., 


the search will end at that final wave number. 


Key Variables Used _in BASELN: 
DIFF (IF) : For the IF-th file, 


(Maximum normalized amplitude value) 


(Minimum normalized amplitude value). 


VALUE : Equal to the input variable VALID. 


(The “validity" factor). 


FE 


IGURE 1 - Example of Baseline Computations. 


Me 
B HT (IP) 
/ Ns verb 


-IPK on ee eS 
/NAVENO (IPK) 


340. 00 320. 00 300. 00 280. 00 260. 00 


TY ET 


Key Parameters Cal_ulated by BASELN: 
By referercirg Figure 1, the following parameters are 


calculated in order to describe each individual spectral peak. 


Variable Name(s) Definition 
1) A and B (Amplitude at the apex of the peak) 
MINUS 
(Amplitude at the minimum point(s) on 
the peak). 
2) IP The IP-th peak in the spectrur, reading 


from left to right. 


3) AB(IP) The baseline amplitude at wave number 
WAVENO (IPK) . 


4) HT(IP) The actual peak height for the IP-th 
peak in the spectrum. 


(Amplitude at the apex) - AB(IP). 

5) IDWORD(I°K) : No peak at data word IPK. 

: Invalid peak at data word IPK. 

: Valid peak at data word IPK. 
Invalid peak at data word IPK, 
measured by the COBASE subroutine. 

: Valid peak at data word IPK, 
measured by the COBASE subroutine. 
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The measurement of peak "validity" is a technique utilized 


to try to eliminate random ‘noise’ within an FTIS spectrum. 
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The method for determining the "validity" of a peak is as 


folljows: 

a) Calculate the values for A and B. 

b) From trese two values, determine which variable has 
the maximum value, then set it equal to C. 

c) Multiply DIFF(IF) by VALUE: set product equal to D22. 

d) If C is greater than D22, then the peak at this wave 


number (WAVENO(IPK))is a valid peak. 
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3.0.4 SUBR» TINE OBR eT 


£3 4% option... subroutine, COBASE will only 
be execute’ when the input vatirsle MULTPK is greater 
than zero and less than seven. By utilizing the same 
general techniques described in subroutine BASELN, 
the logic in this routine will establish a "commun" 
baseline for adjacent spectral peaks. 

For example, to establish the “common” baseline 
for the n-th set of adjacent peaks, the COBASE algorithm 
first searches for NPHOTO(n), the file that contains 
the optimum picture (or example) of these adjoining 
peaks. After calculating this new “common" baseline, 
the data words locating the newly measured peak heights 
in this file also become the reference locations from 
which to measure this "common" baselines' peak heights 
in all the other FTIS files. 

It should be notec that the maximum allowable 
number of peaks per “common" baseline area is twenty- 


five (25). 
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3.0.5 SUBROUTINE PKTABL 


PKTABL is an output subroutine that zenerates two dif- 
ferent types of tabulations listing FTIS peak height informa- 
tion. The first type of tabulation, as illustrated by Figure 
2 , is a "Normalized Peak Height Information" listing that 
is printed for every FTIS file analyzed As illustrated in 
Figure 3 , the second tabulation entitled "Peak Height Tabula- 
tion for all the Files" is strictly a listing, referenced by 
the Data Word location, of all the peak heights found in 
all the FTIS fiies examined. 

A check is also made to determine whether the mechanical 
property correlation subroutine REGRES will be called. If so, 
then all the peak height values that are simultaneously 
located at the same Data Word location and valid in every 
FTIS file examined, will be stored on a temporary disk file 


(File 27) for use as input to the REGRES subroutine. 
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3.0.6 SUBROUTINE PLWAVE 


When requested by the user, this routine plots the 
Wave Number versus Amplitude values for the individual FTIS 
spectral files. Depending upon the value of the input 
variables IPLOT and NMPLOT, the option is available for 
plotting eitner "Non-Normalized" or "Normalized" Ampiitude 


data values or if so requested, both types of spectral plots. 
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3.0.7 SUBROUTINE REGRES 


Subroutine REGRES is also an option that, when requested, 
(input variable MPDISK = i) will determine whether there are 
any correlations between FTIS spectral data and corresponding 


propellant mechanical properties tabulated by the E410 computer 


program. 


The primary logic in REGRES provides a multiple linear 
regression model that defines as its independent variables 
the set of "valid" peak heights written on File 27 by subroutine 
PKTABL. Regarding the dependent variable(s), they are 
determined by the input variable NPHYSP(k), the array that 
specifically defines what mechanical properties the user 


wishes to examine. 


For each dependent variable, a multiple correlation 
coefficient (RMULT) is computed. For any dependent variable 


where RMULT =.756, a non-linear regression attempt is made 


using that single dependent variable and those independent 
variables with non~zero linear regression coefficients. At 
the end of this attempt, a test is made comparing its multiple 
correlation coefficient with the linear one, thereby, 
determining the best “curve fit" for tnis particular dependent 


variable (or mechanical property). 
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4.0 £490 INPUT INSTRUCTIONS 


Input variables for £490 are contained in a single 
NAMELIST called /FILES/. By using NAMELIST, a free form 
input is availiable requiring only that the first column 


on each input card be blank. 


The following table, TABLE 1, presents a complete 
listing and description cf the variables necessary to 


execute the E490 computer program. 
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4.1 JOB CONTROL LANGUAGE FOR CDC SYSTEMS 


The following are examples of the Job Control 
Language needed to execute computer program E490 on 


a CDC-6000 series computer. 


OPTION ONE: RUN REQUESTING NO FLOTS AND NO 
STATISTICAL CORRELATION. 


SEQ, THAOL. 

THAsTLOOSNTL»P4. LAsYOURNAME ,30 

ATT ACHs lL GO, E490G0, IN=YUURNAME, MR=1. 

REQUEST,» TAPE] »NT sHDs Se NORINGs VSN=FTISTAP ENAME 
FILE(TAPEL,BT=K,RT=U,RB=1 gMBL=1257 »MRL=1 257 9 CM=NOy MNR=24,KNB=26) 
LOSET(CFILES=TAPE1) 

LIBRARY,BITSLIB. 

MAP,PART. 

LGO. 

UNLOAD, TAPE. 


Card Column 1 
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Lg 


Bd OPTION TWO: RUN REQUESTING PLOTS 8UT NO 
_ STATISTICAL CORRELATIONS. 
a SEQ, THAOL. 

Ea THA,TLOOeNT1L 994. LAs YGURNAME » 30 

7 ATTACH, LGN, E490G0, ID=YOURNAME »MR=1. 

Le 3 REQUEST, TAPELyNT yHOy Sy NORINGs VSN=FTIS TAP ENAME. 


FILE( TAPFl,BT=KyRT=U,RR=1yMBL=1257 pMRL=1L257_9CM=NO, MNR=245 MNB=24) 
LOSETCFILFS=TAPE1) 


LEPRARY,BITBLIB. 
MAP, PART. 
LGO. 
UNLOAD, TAPEL. 
: REWINL. TAPEB.s 
‘ REQUFST, PLOT sHD»RINGyVSN=PLOT. YOURNAME 
‘? COPYBF, TAPEByPLOT. 
4 RETURN,PLOT. 
WwW 
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REITER SNR mrermenrES! 


OPTION THREE: RUN REQUESTING STATISTICAL 
CORRELATIONS BUT NO PLOTS. 


SEQ, THAOL. 

THA,TLOOgNT1i P42 LAs YOGURNAME, 30 

ATTACH, LGO,E490G0, I D=YOURNAME »MR=1. 

ATTACH, TAPE 10, MECHPROFFILE, 10= YOUR NAME »MR=1. 
REQUEST» TAPEL, iT eHDs Sy NORING, VSN=FTISTAPENAME. 
FILE(TAPEL sBT=KyRT=UyRB=1.MBL=1257_MRL=1 257, CM=NOe MNR=24yMNB=24) 
LOSETCFILES=TAPEL) 

LIBRARY, BITSLIB. 

MAP,PART. 

iGo. 

UNLOAD, TAPEL. 


OPTION FOURS RUN REQUESTING PLOTS AND 
STATISTICAL CORRELATIONS. 


SEQ, THAOL. 

THA, TLIOONT1 P44. LAsYOURNAME, 30 

ATTACH» LGO,E490GN, ID=YOURNAME » MR=L. 

ATTACH, TAP FILO, M-CHPROPFILE, ID=YOURNAMF »MR=1. 
REQUEST» TAPF1,NT sHDy SyNORINGs VSN=FTISTAPENAME. 
FILE ( TAPEL sBT=K,RT=U 8 B=1,MBL=1257,MRL=12579 CM=NCy MNR=24,MNB=24) 
LOSETIFILES=TAPEL) 

LIBRARY, BITSLIB. 

MAP,PART. 

LGO. 

UNLOAD, TAPEL. 

REWINO,TAPES. 

REQUESTyPLOT,sHDsRINGyVSEN=PLOT. YOURNAME 
COPYBF,TAPFS,PLNT. 

RETURN,PLOT. 
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5.0 E490 - SAMPLE CASE 


The card deck shown in Figure 4 is the Job Control Language and 


input data that generated the sample case illustrated in Table 2. 


F49)) SAMPLE CASE 


a a re ne eran ete cae 


SEQ,yTHAOL. 
THASTLOO NTL ye PS. LA SMETHE 29 

ATTACH, LGN, © 49060, 19=SMILTHG»4R=1. 
ATTACH, TAPFLO,PRNOFILETP25, LI=SMITH MR=L. 
REQUEST, TAPEL NT gH 9S 9 NORING s VSN=TAPF25.4 
FILECTAPEL, AT=K »RTH, 2 RaL MAL =1257,4RL=ALIE TZ CMEND, MNP H7G MNR=74) 
LOSETCFILES=TAPE1) 
LIBRARY,BITBLIS. 
MAP,P’RT,. 
iGO. 
UNLOAD, TAPFL. 

$FILES 

ITAPE="TP?5', 

ITOTAL=4, 

INFILE=16920,74, 28 y 

MULTPK=1, 

NPHOTO=24, 

LLIMIT=3030., 

RLIMITHC TIN, 

WAVNOR=2 3500, 

VALID = 2.92, 

MPNTSK=1y, 

MPPECS=17,2142552S, 

NPHYSP= lee, 3e4y 


TPLOT = Oy 
NMPLOT = Oy 
$F ND 


FIGURE 4. E490 SAMPLE CASE 
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6.0 E410 ~ COMPUTER PROGRAM DESCRIPTION 


£410 is designed to create and maintain a permanent 
file for solid propellant mechanical properties. Once 
established, this file will be utilized for optional input 


to the E490 computer program. 


This FORTRAN code consists of two parts, a main anda 
subroutine called RENAME. The main program creates the initial 
file and prints out data contained in the file. RENAME is 
designed strictly for updating or accessing an existing 


propellant file. 
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7.0 F410 - INPUT INSTRUCTIONS 


E410 is always executed utiliziig one of two possible 
run modes; the CREATE mode for building a new propellant file 
or the UPDATE mode designed to update or access an existing 


file. 
7.0.1 CREATE MODE 


The first card in the input deck must have CREATE in the 
first six columns. This tells the computer to stay in the 


Main program. 


The second card is the title of the chart that is printed 
when the data set is established. The title can be 40 alpha- 


numeric characters long. It is read in on a 10A4 format. 


The next sequence of three cards is repeated for each 
set of data. The first 15 columns on data card type 1 are 
ignored. Any comments with alphanumeric data can be placed 
here. The next five columns are ditadated to the specimens 
age and the next two sets of ten columns are allocated for 


the aging temperature and test temperature respectively. This 


card is read with an (A8,7X,1I5-,21I119) format. 


Data cards 2 & 3 are used to input the values for modulus, 
strain at break, strain at maximum stress, maximum stress, strain 
energy density and strain endurance respectively. Ten columns 
are allocated for each of these plus four extra variables that 


can be allocated later. Five variables are read from each of 


these cards. Only the first 50 columns are read on each. Ail 
five fields are end to end. The format is (5(F10.3)) for both 


cards. 


7.0.2 UPDATE MODE 

In this mode, the first card in the input deck must have 
UPDATE in the first six columns. Also, tht title card is 
eliminated and the first of the three data cards 1s changed. 
Data card 1 in this mode should contain the action to be 
performed in the first 8 columns, skip the next two columns, 
and have the line number that the action is to be performed 
on in the next five columns. These two new inputs replace the 
‘columns ignored on the first data card in the CREATE mode. 
The remainder of this card and the two following data cards 


are the same as in the CREATE mode. 


All actions should be nested together by type, and these 
groups then input ONLY in the following sequence; 

1. All CHANGE cards. 

2. All ADD cards. 

3 All DELETE cards. 
If there are ro CHANGE cards in the input, then the proper 
input sequence should be all the ADD cards first, then any 
DELETE cards next. This correct sequencing of the input 
cards is mandatory in order to maintain the integrity of the 
existing properties file. Care shou? also be taken to put 
ihe proper data cards after the CHANGE and ADD action cards; 
the other acticu cards do not require any data other than the 


action and line number. All strings of actions need to be 


followed by either ENDLIST or ENDNOLST. 
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8.0 E410 - SAMPLE CASE 
mt EE, CAE 


The card deck shown in Figure 5 is the CDC Job Coitrol Language 


and input data that was used to create the propellant properties file 


listed in Table 4. 


Card Column 1 “410 SAMPLE CASE 


SEQ,THAOL. 
THASTLOO2P4. LAVSMITHO, 39 
ATTACH sLGNsF41 060, IN=SMITHO,¥R=1, 
MAP,PART, 
REQUEST, TAPES, xoF, 
LGO. 
CATALOG, TAPER, TOHE A 56V612,1D=SMTTHORP=990, 
CREATE 
TEST CASF --~ SEPT, 27, 197¢, 
ADD 4 77 77 
1119.0 0.0 24.00 165.4 0.6 
0.9 2.0 0.0 0.0 0.0 
ADD ? 190 77 
0843.9 9.0 21.96 126.7 0.0 
oO." 9.0 920 0.90 0.0 
ADD 2 170 77 
774.0 a) 22.46 171.2 0.0 
92.9 20 9.9 0.0 92.0 
; AND 2 150 77 
B45.0 9.00 71.090 116.1 De 
3 0.0 909 9.9 0.0 0.4% 
ADD 2 139 7? 
= 773.0 0.9 24.3 117.8 0.0 
i 0.9 a) 9. 0.0 0.0 
E Ann ? 110 77 
| 932.9 9.0 24.3 144.9 0.0 
f 9.0 0.0 0.C 0.0 9.0 
: ann ? 75 7 
i 1171.9 9 Le 151.8 0.0 
k 0.9 9.0 0.9 6.9 6.20 
t ANN 4 1930 77 
E 142, 9.0 15.8 156.2 0.0 
b 0.0 0.0 0.0 0.0 0.0 
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N9ANn297¢ 
IOAN 29R{ 
99097 9C/ 
ONOd AE! 
INOVIDIE 
9000392¢ 
900030231 
0990394: 
0007305! 
0000306: 
900030 TE 
GONNA VWAC 
FONG209( 


OBONALL: 
non0212! 
700031 3¢ 
00003146 
nonaass! 
NOODALE! 
9990317 
99003'1" 
N9N0319¢ 
900037* ¢ 
DONNIE 
BING +724 
9000 42 2¢ 
5000324( 
06N0%25( 
300037€! 
N1000%7 7 
000032AC 
0600329 
0000230¢ 
90002717 
00002 52¢ 
90n03 aar 
N0N02346 
SN00335f 
0000374H¢ 


Pe Ae FAMERARSKER ESC SAKREEHE SEKERERSHESEE CHS ES EE EKEE SHHESCSHEHHEEEHEENIO NI 2 TC 
of HECK HOW MANY PHYSTCAL PROPERTIES ARE BEING EXAMINED. 


NONOIEAL 


0000349 


c 
HUMANE - 5 0060356 
NUMIND = COUNT 9000351 
NN 625 13=1,10 0000357 
TE CNPHYSP (13).NE.0) NUMDEP = NUMDEP + 1 0000353 
NPHYSIOT3) = NPHYSP(1 3) 00003%¢ 

525 CONT INGE 6090356 
NUMOBS = ITOUAL 0000355 
Nn 630 LL= i, NUMOBS 0000357 
MPFILFCLL) = MPRECS(LL? 0000358 
630 CONTINUE 0000359 

' 9000360 

Ceocecwoeol Att FIF MUETIPLE LYINFAR REGRESSIGN ROUTINE: 9009341 

C 0000362 
CAGE Br GRES 9000363 

¢ 0000364 

Ce eRAS HAAR ES KEAEEKHE TEKH SEEAEER SESE CESS EHOSAH ES. CSE y ESCSHSH OS CESS HS SSS 090365 

¢ NONOIKA 

650 TE CLOLITLFQ. LOR. NMPLOTLEQ.L) CALL PLOT (00,0.59999 

C 2999369 

{ ,escenco el ORMAT STATEMENTS. 90003 75 

c o000nsTi 

TIS FF IRMATUA2» Ag 3A27 — Ke SAZS 0060372 
T7270 FrOMATISSX.*NUMBER NF COMMIN BASELINE AREAS SPECIFIED 0.*) 0000373 
7725 FORMATI/5*."NUMAFR OF COMMON AASFLINE AREAS SPECEFFED ",115%.70900374 

i} 0009375, 


73S 7 VRMATELOK,FAPEA M641," TS BEST PICTURED IN FILE "9TBy% 2% p%e 1 0X,” 
TEE GEMTT ESTEMATE = "sFT.2e/,LOX,PRIGHT LEMEY ESTIMATE = ", F722) 

P25 CARMATE SSX, OND STATISTICAL CORRELATION WITH PHYSICAL PROPERTIES W19690378, 
TEL PE PERE ORAENS C4" DESK=000%) 900379 

Ya 7 IPMATO/SKS"STATESTICAL CORRELATION WITH PHYSICAL PROPERTIES HAS 80009380 
ITTN REQUESTER "//5Xe"PHYSICAL PROPERTY INPUT DISK  FT®,12,"F00L."0090381 
») 


} 99N9382 

f THeMAVESSXS* PHYSICAL PROP, RECORD NUMBERS 7, 1001 3¢S5K9736X,_, 10083, 80999283 
rs FOX LOOCTI,SRD) NONNIKRSG 
eR CAV ATE “SX, *® THE FOLUOWENG, PHYSTCAL PROPERTZES WILL BE USED™/5K, 7450090385 
TOF NOP NYT VARETABLES ®) 6009336 
TA SC URRMATEIGK,*PROPL NO. ToTE Ss? — @,6A4) 0000387 
7617 EC TIMATE LUE, *%, @SUMMARY OF AMPLITUDE NORMALT7ATION "// 9900388 
Loreto SP PETL YM TSK, ME EL EMSS 5K, 3( ONO, EBX) - PND) 0000389 
toe “WRGATISGX, FEL BeLTXP IZ) 9000390 
RAP FHRMATI S/T. PMAKT MUM AMPLITUDE NEAREST *,F 6.19% We = 9, &FFL5.7,5900039) 
1X97) 900039? 
cw? FIQMATS TU," dAVE NUMBER AT MAX, AMPLITUDE = (WNMAXS = %,F10.4, 30100900393 
tXeFENLSIS) 90003946 

"2% FVORMATECOX,PRASELINE AMPLITUDE AT CWNMAX) = 9,8X,4CEU5. 725K) 7) 5900395 
RAS FORMATE TX, @PFAK HEIGHT AT (WNMAX) = oe 1 5X,60F15.7 9 5RSA) 3099396 
RAN ENAMATETL, “CLXBE LOD) 0990397 
FOO CNAMASELXs FOSKeFLOL4)) 0000398 


FAD FORMATIOX GS 134K oF LO. Ge PKeE LS Tc 20K es F359 4K eo FLO. Se PK ELSA TDD 7090399 


940 EORMATCLHE s4Xe"FTIS FILF NUMBER “,13,20X,*MON-WORMALIZED (PURE} $00004 


UPFCTRAL DATA™//5X,_ 30 "DATA@ 8X9 @WAVE™ 2, BK» PAMPLE TUNES 8K) / 90004 
75X— TOMHIR Gs TX, "NUMBER®, 28X07) 00004 
74t FORMATCULHL 4X, "FTTS FILE NUMBER °,93,20X, "NORMALIZED SPECTRAL DATO0004 
LA" //5X_ FCP DATA BX, “WAVE™ 8X, PAMPLITUDE@ ,OXS/ 00004 
25K, ZEPHORDM, TX, *NUMBER™ 24X93 739 00004 


950 FORMAT(T<s"NORMALTZEUNG FACTOR FOR THES FILE = "<SXe4(EL5.7,5X)) 00004 
F7N FORMATCINL SX,"FNURTER TRANSFORM INFRARED SPFCTROSCOPY COMPUTER PRO0104 
LOGRAM™/5X,*E490 — THIDKOL CORP. / MUMYSVILLE . ALABAMA 35890707 7//5xK00004 


2e* INPUT DESCRIPTION "/7) 00004 
972 FORMATC/SX "NAME OF FYIS INPUT TAPE "Aa? 00004 
375 FORMATC/SX,"TOTAL NUMBER OF FILES CTHIS RUND *,53-%.*) 00004 
DRO FORMATE/5X e*LIST NE FILES TO BE ANALYZED *, LOC13, 5X) 736X, L001 3, 500004 

1X9 /36X,_ LOCT3,5X)) 00604 
JRL FORMATC/5X, "NORMAL EZING WAVE NUMRER #F6.1, M8) 00004 
9B? FORMAT(/5X_,"*NON-NORMALIZED SPECTRAL PLOY HAS BEEN REQUESTED. *) 00004 
IRI FORMAT(/SXKs®NORMALIZFD SPFCTRAL PLOT HAS BEEN REQUESTED. *) 00054 
78S FORMAT{(/5X,°ND SPECTRAL PLOTS HAVE SEFN REQUESTED. *3 00004 
995 FORMAT(/5X,"°F AK HEIGHT VALIDITY =",F4.1,"% PER CENT."} 00900 

sTnP Oo 00004 
ENN 00004 
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MWaeOU eee APCALC CAMPE 10, KAVENG ,WAVNOK, .~ MN, JAM, APHEX WE, N12, AP2) 


arr a SEPTEMBER &, 1977 - CDC REVISION. 
C 
eee eUOGaC WELL CALCULATE THE NORMALIZACEON AMPLATUGE VALUE CAPNORK) 


Cea... =... FOR THE PEAK NEAREST WAVE NUMBER (WAVNOR) CM-2. 
C 
COON SSPECTR A LTETLEC30) ~NPT29 €30) 
CHMMON /TMENPK/ ALEMITC6),ZLIMSIT(6)_-SFILE(6) .NORACO.NUMGEr 
ULMENS:ON AMPLTN (468) ,WAVENO(468),A025) 
t=0 
}-9 
NU 1 s328,25 
{ A(@J3b = OC. GED 


C weaves SE4RH FIR PAXIMUM AMPLITUDE NEAREST TO Wi¥NOR. 
5S fe Toe 

TE (WAY ONOCT DCT. CWAVNOK420.0)) GO TG 5 
tr TWAVENC(L).LT. CWAVNOR-20.033 GO TO 1G 
it €5,FQ.0) Ti=a -} 
be 343 
ACS) = AAPLIDOTD 
5) Tt 5 


pm AP YAY AY) 
ms 1S Jbe7,d 
it (AOCTITLT Oo OMe Xb GO TO 15 
APMAX = Av JJ) 
TRFF mv) 
‘7 ONT ENUE 
MK = PEterRes 
OPT SVOTEE = ik 


(.+ 'OK st See Pf A “COMMON BASERIP! §S Tu BE USED 
Ce AR ie + THE MORMALTZAcCEON AMPLITUDE.. 


it (NOON EG #1 GO Tl 45 


ee Bh oC are MINTMUM® POINTS ON BOTH SIDES OF APHAX PEAK. 


x7 M+] 
7 UF fader TDRKDL CY AMPLIDIK*1)) GO FG 340 
a= J 
5Uu TA 290 
3 RMIN = AMPLIOCKELD 
KMEN = K41 
: 
(Non 21 


+S FE CARPLID(K+1L).GT.AMPLTD(K)) GO TO 40 


ee ae ee A RE OF 


K = Kel 
GO IN 35 
40 RMAX = AMPLTO(K) 

KMAX = X 
GO TO 100 

C 

Leccecceel = THE I-TH DATA WORD IN THE SPECTRUM FILE. 

c 

45 1=0 
200 I=F+l 

TF C(WAVENOCIT).GTALIMIT(JAXI}GO TO 200 
is = I 


710 TF CAMPLTOCI-159.GTAMPLTDO(I}) GO TO 220 


BMIN = AMPLTD{I-1L) 
KMIN = [-1 
Y=I-1 
60 TO 210 

220 f=IS 

230 IF (WAVENDUT 41) .LEZLIMIT( JAX} I60 TO 240 
f={+] 
60 TO 23¢ 

740 f=] 

245 TF CAMPLTOCT+ID GTAMPLTO(ID) GO TO 250 
tr Pel 
60 TN 245 

250 RMAX = AMPLCTOCES 
«MAX = { 

C 
iftr CHGX = WAVENOTKMAX) ~ WAVENOCKMIND 


CHGY = BMAX & TN 
£ 
Cav sccee CHECK FOR CHGY = 0.0. 
{ 
TF (CHGY.EQ.G. EO) GO TO 110 


C 
Cr ewcenne MULTIPLE TANGENTIAL CHECK. 
c 
{fF (BMIN.GT. BMAX) GO TO 80 
t, 
Lie ec cane RIGHT = SHORT SIDE OF THE PEAKES). 


75 SLOPEL = €(CHGY/CHGX) 
CY = AMAX — AMPLTD(KMING1) 
CK = WAVENTCKMAX) — WAVENO(CKMING1D 
ScOPt2 = (CY/CX) 
SIGN = SLOPEL * SLOPE2 
ZZ ABSCSLOPEL) 
ZZ2 ABSCSLOPE2) 
BMIN = AMPLTD(KMINGL ) 
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CHGX = Cx 

Coy = CY 

KMIN = KMIN ¢ 1 
60 10 75 


v 
CoscuceeeLFFY = SHORT SIDE OF THE PEAK(IS}). 
€ 
40 SLOPEL - €CHGY/CHGX) 
CY = AMPLTOUKMAX-L? — BMIEN 
CX = WAVENO(KMAX-2} ~— WAVENOCKMIN) 
SLOPF2 = &CY/CX) 
SIGN = SLOPE] * SLOPE2 
ZZa = ABSESLOPEL) 
£Z2B = ASS SLOPE2) 
TF (SIGNLT.O.0EO-OR-ZZBLT2ZZA} GO TC 130 
BMAX = AMPL TOC KMAX-1) 


THGK = UX . 
CHGY = CY F 
KMAX = KMAX -— li 
60 19 20 
L1G DI = (CHGY®( WAVENO (MKD -WAVENOCKMIND DD SCHGX 
2 = DL + BMIN 


AP? = APMAKX - D2 
NE = dAVENOCMK) 


en eeeaS END THE VALUCS FOR (APMAX,WE,D2,422) BACK TO @AIN. 


f 


RETURN 
enn 


I 
| 
: 


we ee etna, 


SUBMOUTING BASELNOOLEF, VALUED oOnndGgooF 


" 9000007 
L - oe SEP TEMBIR 1391979 ~ REVISION FOR COC 6600. 000G006 
t 9000009 
CUMMON /POINTS/ NWORDSU30} 0000010 
NIMFNSION AMPLTO1468) yNAVENO( 4668), IDWORD(468 } NoONOLL 
IMENSEGON CHINI (50) ,CMIM2050) pPEAK{ 50) 2ABC 50), 47( 50) pDIFE(3C) 0009012 

LOK ARERR SHOE SEER EEEREEETEEEERE SSF HSESS SOSOEES € ETREVYIR OES ED EEE OE? EO 00001 3 
rq o000 014 
Core vent lS SUBROUTINE SEARCHES FOR AND RECORDS VALTO PEAKS 0000015 
Cian -weeWI THIN FACH INDIVINUAL SPECTRUM. ooooole 
3 00000! 7 
PO AR RERKMT OH AE EE TERE REEEEEEEEEKER EE SEE SEES REESE EE HEE REAERED SOS HHS REL EDONNNTS 
C 0000019 
Casccas efltcte FOR MAX. 50 PEAKS PER INDIVIDUAL SPEC TRUM, 0000020 
C 0000021 
Uae ~. .f = tHe I-(1t DAtA WORD IN THE SPECTRUM Fiitc. 0000022 
c 0000023 
an » et? - TKE YP -1H YALIN PEAK TN THE SPECTRUM €1 EFT-RIGHT). 00000°4 
t 0000025 
Pas cee ee POWORNETOKD = 0 NO PEAK AY DATA WORD ""’PK*. 9000026 
t 1 INVALIO PEAK AT DATA WORD *IPK". 0000027 
r 2 VALID PEAK AT DATA WORD ®1PK®, 0009028 
\ 3 INVALID PEAK AT DATA WORD ®7PK" , HEASURED0000N29 
c BY THE COBASE SUBROUTINE. 0000030 
é, 4 VALID PEAK AY DATA WORD FITPRK® , MEASURED 000003: 
f BY THE COCBASE SUOSROUTINE. 9000032 
% 0000033 
CERES RECEE TAHA E ERNE EEETESEREREESRS FE ORESSSSSSS SORSRERG OE SNS HHS SEF = 0000034 
5 0000035 
See, need EUTALi2* SOME VARLT ABLES. 0000036 
7 9900037 
rare 0000038 

sa 0000039 

wes 0000040 

WY DO + - Ly GbR 0000041 
POWHKA Lh = O 0000042 

er COAT ERE 0000043 

WX TION Fs {8C 0000044 

Ast { i ’ Oe U 0000045 

eu¥d fb = GLO 0000046 

#tor CONTINUE 0000047 
t=O 0000048 

€ 0000049 
COE Ree ER EEK ERERE EES + EEEETERESEESES RS BESERSES HD $F SHOPS SHVSES EOE SSS 46460000050 
f, 0000051 
REARIZ1) 1, AMPLTD,WAVYEND 0000052 

NW - NWORNSTIF) 0000053 
0000054 


“G Tot 0000055 


oe a 2s ae 


lt (WAVENOUL)GT.3700.E60) GO TO 10 

15 TF CWAVENOCT).LT.2700.£0) GO TO 90 
it (4.NE.09 60 TO 20 
CMENLCIP) = AMPLTOCT) 
Nl=1 

20 TE (WAVENOCL EL) LT. T000EO.0R.(101).GE.NW) GO TN 130 
1F CAMPLTOCYE+1L)2.GT.AMPLTOCED) GO TO 30 
ff (5.NE.0)3 GN TO 40 
CMENLCTP) = AMPLYOCI4L) 
Ni= T#l 
T= T+! 
Gf Ta 20 

30 J=1 
f=Uel 
Go a 29 

40 TOWEN=! 
WRETE(6~400) LOWD Nea MPLYOCI) sCMIRLGIP) 

400 FORMATOIX, P55 3XeF 15. 7e3KvE15073 
A = AMPLTONCL) ~— CMINLCIP) 
PEAKCIP) = AMPLING 3 
Pn = 1 

50 Tt CAMPLIOLE +1) GT. AMPLIDEID) GO TO 60 
V=Fel 
sn to 50 

£0 CMIN?CIP)= AMPLTD(T) 
NZ=1 
MIN? = N2 
BMP2 - CHIN2ZCTP) 
\ - PEAKIEP) - AMP? 
roo: AMAMICA,S:? 
f.¢2 - DIFFCYS) © VALUE 
We €C.GT.922) TOWORDCTPK) = 2 


4 3OVAL TN PEAn CQELOW eeee 
CALCJLATE FQUATLON FOR LINE CONNECTING TWO MINIMUM POINTS. 


ANJUS. BASFLINF FOR THD YANGENTIAL CONDITION: 
THUS BASELINE MEINENUSN POINTS WILE LIE ON THE TANGENT. 


racy 
«tay 


AAVENCENZ)D - WAVENOCUNL) 
CMEN2Z(1TP) ~- CMINECTIP) 


uw 


CHECK FOR A = B= C. 


{F (CHGY~EQ.0.0F0} GO TO 110 
TF (C.EQ.8)9G0 TO 80 


3 = SHORT STOF OF THE TP-TH PEAK. 


7> SLOPE! = CCHGY/JCHGX) 


0000056 
0000057) 
000005a' 
000005% 
0000000! 
0000061 
0000062: 
0000063! 
0000064: 
0000065: 
0000066! 
0000067! 
00000681 
000006% 
00000701 
0000071( 
00000721 
000007 3¢ 
00000 74< 
oo00075« 
0000074(< 
0000077( 
0000078 
0000079( 
0000080( 
000008) ( 
000008 2( 
6000083( 
0000084( 
9000085 
0000086( 
0000087 
000008aC 
0000089¢ 
00000905 
00000911 
0000092 
0000093¢ 
0000094: 
0000095¢ 
0000096t 
0000097 
0000098 
0000099 
0000100¢ 
0000103 ( 
0000102¢ 
00001036 
0000104 
0000105¢ 


F 
F 
i 
f 
t 


Pre Wey ore 762 FPO TP 


wh 


re OE TT GPR, TEE ER oro > 


aa wet 


aKeeigoansss 


: 
E 1 


90 


Ct - CMENZGTP) -— AMPLTDEND4LS 


OO00LO46\ 


CX = WAVENNTON2) - WAVENO(NLSL) O0000107C 
SLOPE2 = ICY/CX) 0000108¢ 
SIGN = SLOPEL * SLOPE2 Qo0oLost 
771 = ABSESLOPFL) 00001106 
722 = ABSCUSLOPE2) OOOGLIIC 
LE CSIGN.LT.O.0£F0.ORL7ZZ2.LT.2ZZ1) GE TO 110 o000lr12¢ 
CMINILOCI®) = AMPLTCCUNIT4L) OOOOLL3E 
CHGX = CX 0000114( 
CuGy = CY reLeleles @ axe 
Ni = NE ¢ 1 COOOLI EH. 
Cy Ty 7S oOoOoLLT. 
OGOO1L1&« 

A = SHORT SENE OF THE IP-TH PEAK. ooo0oLl: 9 
00001201 

SLOPEL = €CHGY/CHOX) 0000121« 
CY = AMPLTNEN2-£) — CAMINIT CIP) 0000122( 
CX = WAVENOOUN2-L) - WAVENOCNLD 00001231 
S19PE?2 = OCY/EX) 0000124 
SIGN = SLOPEL ® StLre? 00001 25¢ 
771i = ABRSCSEGPEL) 6000126: 
272 = ABSTSLO": 2} ooool2T< 
Ir (SUGM.t T0290. Re Z7Z20LTO7ZL) GO TO 110 00001281 
CMENZETE) = AMPLTNCH2-L) 0000129% 
CHuGKX = CX 00001 30( 
cacy - CY 00001311 
NN. 7 NS ~ 2 0000132! 
GU Tt 00001 33¢ 
0000134! 

(Civ vee AVEN FS PKE-WAVENOCQNE DF) SCHGX 9000135! 
00001 36( 

ARCLLP) = AMPLITUDE OF THE BASELINE FOR YHE IP TH PEAK. 0O00L 37! 
HT( P30) = HETCHT OF THE IP-TH PEAK. 00001381 
00001391 

2 ES p> OD ECMINE CTP} 00001401 
areca CEAKITP) - ABCIP) 00001411 
> Cy Dw VC EPK)D.FQ.2) GOTO 120 0000142! 
re WNROCR Yt, = 3 0000143) 
ivi ype} OOOO 44: 
CO0OL4S: 

(AMP 2,EMINZ? THE TRUE MINIMUM POINT ON THE RIGHT SIDE. 0000146: 
000014: 

i @ENTCEP? = AMP? 000014581 
NI=UMIN? 0900001491 
t=1+4) o0001Ls0 
GA TN tS oo00151Li 
0000152" 

tk (WAVENN(1).LF.1800.£0) GO TO 100 0000153! 
I= el 0000154! 
oy ™ 99 0000155: 


111 


100 TE (WAVENDCEDL1. 700-E0.0R.1-GE-NW) GO TO 130 0000156( 


IF (K.NE.0) GN TO 20 0000157¢ 

J=0 0000158 
CMINLCIP) = AMPLTICI) 0000159¢ 

Ni=I 00001606 

Kz] 0000161¢ 

co 9 20 0000162{ 

1B 00001630 
CEReRERERAREAEKEREEEREERS ONELEEEREEES SEERESEE SO SSS EE OL EEE RELE FES SSEES SEES 1000 164( 
c 0000165¢ 
© vecnee WRITE FTOWNRD TNFORMATION TO THE DISK FILE 22. 00001665 
C oo0ooL16rec 
CEGREEED EE RERE KR CEREKIZE SESE SSESEESSESSHVKS CES SKKSSESSL SECS HESEHES OTE SD S44 ¢ 90001 691 
C 0000169¢ 
130 WRLTE(S, 5000) TF 00001 70¢ 
WETTES&,5091) FOWORD 00001 71L¢ 

5009 FLIRMAT (5K _15) 00001 72¢ 
SOOL FORMAT €50/5K, LOOT LD) 0000173¢ 
ARITE (22) TF, INWORD 00001 74C 

WRITE (24) TF ABH: 0000175¢ 

RE TURN 00001 76C 

END 0000177TC 
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SHUAKOULTINE COBASE COTFF » VALUE 5 JX) O0DNG: 


C 00000: 
Cec ew se e SEPTEMBER 28, 19796 (CHANGED AT AF/®PL =: QCT. 24, 1979). 
Cae weewee 000001 
CoscwecesfLDC REVISION. 00000) 
C 00000! 
COMMON /SPECTR/ TF ELE€30) oNPT29(30) 00000) 
COMMON /STWINPK/ ALIMIT(6),ZLIMIT(6) »MFILEC6) »NDRMCO, NUMBER 00000) 
NEMENSTON AMPLTD(468),HAVENOL 468) 5 IDWORD( 468), 28650) ,HT E50) 00000) 
DIMENSTON EDGEI( 25)», EDGE2 (25) _1 SPEC (29) » JNORD( 25) »RASTERE 253 00000} 
NIEMENSTION A( 2599 Y€50) ¢2€50) ,OEFFE30) 00000! 
re ee ete ee ee ae eee ee ne ee ee in en er = 900001 
Cc 00000) 
Ceosrceese efkOMMON BASELINE ROUTINE / MARCH 19, 1979. (ocs) 000001 
Cc 000002 
Cue-csceeTHts SUBROUTINE IS A USER REQUESTED OPTION. 00000. 
Con-ceceseTHIS LOGIC IS uIILIZED TO ESTABLISH A COMON BASELINE 000007 
CoeseacvcFOR ADJACENT SPECTRAL PEAKS. 000002 
Cc 000902 
{an -- + ---- ~~ -- a ee ee 3 000002 
C 000007 
Crs veweeel - THE [-TH DATA WORD IN THE SPECTRUM FILE. 000032 
C 000002 
Coeeceeee IP? = THE 7)-TH PEAK IN THE COMMON BASELINE (LEPI~-TO-RIGHT). 000002 
c 000003 
Cen sese. -MARIMUM NUMBER OF PEAKS PER BASELINE = 25. 000003 
c 000003 
Cosas oe eAY MUM NUMBER OF COMMON BASELINES PER COMPUTER RUN = 6. 000003 
c 000003 
Co see ee se ee ee en a ~~ - - 900003 
c 906903 
DQ 2 222,29 000003 
USPECLI2) = 0 000003 
? CANT INUF 000003 
c 000004 
Cen esee cs efECIN BY SEARCHING FOR THE FILE (MFTCECEX)) WETH THE OPTIMUM 900004 
Cece anee ePIC TURE OF THE MULTIPLE PEAKS. 000004 
f, 090004 
B=) 00000 + 
IN 5 £ - geNUMBER 00000+ 
ie CLFILECT).FQ.MFILELIXD) GO TG 4 000004 
ISPEC(13) = f 000004 
iA = th #1 009004 
60 1 5 000004 
4 1F = 9 006005 
5 CONTENUE GIGONS 
MAXX=NUMBER - 1 000005 
Lane = 1 000005 
c 000005 


C ECAD FROM THE NORMALIZED FILE 21. 000005 


| 
b 
| 
i 
F 


SORTA FEE TIE FTE RRR ITEP LPT EE UTERUS TO 2M ANT TN SPE 


Cc 
6 READ (21) TX, AMPL TD, WAVENO 00000 
READ £22) EX,IOWORD 00000 
READ (24) 1[X,AB,HT 00000 
if CIX.€Q.1F) GOTO 7 rerererezey 
C 00000 
C. WRITE TO UPDATED TOKORD FILE 236 00000 
C oo00nno 
WRETE £23) UX, fDWORD 00000 
WRITE €25)3 IX,AB,HT 0000% 
GO TN & 00900 
Cc 00000 
7 1=0 00000: 
tPxi 00000 
VALE = 0 00000 
10 T=] 41 00000 
IF iWAVENO(2).GTALIMIT(05K)} 60 TO 10 00000 
1s = [{ 00000 
20 1F CAMPLIDOEI-L).GT-AMPLTD(I)) GO YO 30 00000 
FOGEL(C TP) = AMPLTOCI-L) 00000 
N= f-t oo00c 
l=t-1 00000 
on TA 20 00000 
9 T=TS 00000 
45 TF (AMPLTOCECLI.LT.AMPLTD(ID) GO TO 40 Go000 
t-te 00000 
LO TO 35 00000 
40 IWDROCIP)Y = I 00000 
ALFET = AMPLITD(I) ~ EDGELCIP) 00000 
55 YF CAMPLTOCT#1L).3TAMPLID(I)) GO TG 60 000006 
2-F+l 00000 
Gte TO $5 00000 
Cc Q0000 
Ce. . pw een TAKE INTO ANCOUNT THE POSSIBILITY OF A DISAPPEARING 2-ND PEAK. 00000 
Coecnsrea WULY « 1979. 00006 
e 00000 
60 fF (WAVENOTE)LEZLIAMIT£{IX}} GO TO 600 oo00oc 
EDGE 2CIPY= AMPLTOCI) 90000 
c 00000 
Ceaec ere oF INISHED DEFINING THE FIRST PEAK. 00000 
o 00000 
62 1P<IP4l 00000 
COGELCLIP)=ENDGEZTTP-1) 900000: 
C 00000 
65 1F CAMPLTIDCI41L).LT. AMPLTD(!)}) GO TO 70 00001 
T= 31 reer ep § 
GO 10 65 00001 
TO SWORDCEIP) = I ocool 
TL TE (AMPLTDCL41) GT AMPLIDUI}) 60 TG 75 o000! 
cooor 


f=itl 


Q900U 


| 

E 

; 

: 1 Yas eleleree of) 
4 7> LUCE ET) = ABPLILGA5 000016 
3 y2 2:4 oocol: 
(F CwsicMias Pobbwce this (okb) GO iG fe 900010 
GO TQ &2 000011 
q t o000e0l! 
G0 20CT2ENPELK? © AMPLIDEDD 00001! 
W2 - 45 09001) 
rp = NEES 000011 
12 = HASTERVIP 000011 
Vin 65 000011 
ooool! 
C ose ee PAPC ATE THE EL EMGTH OF THE RIGHT SIDE. 000011 
j t 000011 
tin ad SMD aD 000012 
Bo bivdd AAPL -OdLZ) - EDGEZCTPD 000012 
USN = fin A LCAL EEL A HED 000012 
| no 7 OPFECAF) * WALUE 000012 
te «H.¥.96.000) GA 19 LOO v0001L< 
‘ 000012 
| ea ae 000012 
u oo00! 2 
Gene) eC ALO TAR £3 TTYOM FOR LINE CONNECTING TWO MEINT@UM POINTS. 000012 
c 000012 
fee eAN y DASEL ENE FOR THE TANGENTIAL CONDITION 3 000013 
Cee ce oe SO BACOLING MENT MUM POINTS WELL LIE ON THE TANGOQN/. C0001: 
: 000013 
: OolNG~ AY ENGI‘ ¢D -WAVERGOND) 00001: 
: ‘HEY -£9G°2CTP) - COGERCL) 000017 
; ‘ 00co!? 
4 wo EHEC OE CHEY = 4.06 000013 
: 00001? 
PF LCHCY £45 OFO? -A TO LLO 00001: 
i 00001: 
i. PEEPLES TANGO YTUAL CHECK. 000014 
e 000014 
iF dou" 2440.52. euGe2trea) GA TH 120 ooo0l 
000014 
; eee tt “.tu.) SELE OF THE PEAKIS). o0col4 
000034 
: LP NE GPRE > CCHGWACUGKD 000914 
CY = FLGE2ZCEP) - AMPLTOCANL#2L) 000014 
X= WAVENOEN2) -» WAVENOENE+L) 000014 
SLaPe? = (CV/CX) 000014 
: SIGN = SLOPEL * SLOPE2 0000;¢ 
{/Z1 = ABSUSLOPEL} 000015 
#72 = ABSISLNFE2) ooodis 
rE (SIGNALT.0.0FO.CR.Z22-LT22Z1) GO TU 110 00001* 
tUGELELD = AMPLIDONLEOD) 000015 
¥ CHOK = CX ooogls 
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CHOY - CY 
Ni. ‘44 4 L 
Gu FA Te, 

& 


Cowcwes se cALEET ~ SHORT STVE GF THE PEAK(S). 


c 
L7G SE NPE, “CHOY CHER) 
CY ANP: @:ONI-2D) - ECGELCL? 
CX = WAVEROOINZ-T) - WAVENOC(NL) 
Stnpe2 = (CY¥/CX) 
sO = StOPEL * SLOPE? 
Jit + tRStSLOPEL) 
442 = ABSESLOPE2) 
{fF @SIGN.L T.C.0EO,.OR.ZZ2eL 7.221) GA Ta ikO 
CUSEZCIP) = AMPL TOI N2—1) 
LiGX = 16 
CHOY = CY 
N2 = N2 ~ 1} 
50 tT ten 
Wty re ttel eae lEPO-OUE PRe OUS PEAK REIGHT INFORMATION “HAT WAS 


Cocres ec eCALCIN ATED IN THE @BASELN® SUBROUTINE. 
c 
Come e ee oS > NG. UF PEAKS FOUND IN BASELINE INTERVAL BETWEEN 
(eee cewiee NL + N2 BY BASELN SUBROUTINE. ; 
C 
169 (SUM og 
Cee REAL) CO Te Lad 
VASNE =D) 
MA SN? =N7 
fe. vo LOS LA. HASNLEMASA- 
IF CLIDWORDILAAD.NE.0) TSUM = ISUM ¢ i 
TOwWORD(LAX) = 6 
1°58 + ONTENUE 
fF CPP WELL) GO TG £15 


DU Lid Ate-~L,lP 
MASTERQNZ} = SWORE SZ} 
CONTINUE 


NPEAK = EP 
RLS PSHITFTENPEAK- TS 


1m 150 J=t »NPEAK 
{Air = MASTER( J) 
NEN (CHG Y* (WAVENO(140)- WAVENDUNL DD} /CHGX 


HL 


eq 


Case ceee cAB( IP) = AMPLITUDE OF THE BASELINE FOR THE 19TH PEAK. 
¢ 


AftJj’ = O03] + EOGELEL 


oOvooL 
ooonl' 
ooool! 
ooool' 
oooole 
0000: 
o0coLe 
oo0ols 
0000 lt 
ooo0ol: 
oonale 
oo0ch 
0000i¢ 
oooole 
00001" 
00001; 
00001; 
00001 
00001" 
000011 
00001; 
00001; 
00001; 
00001; 
o000tt 
oooOolt 
N000Lt 


OooOolE 
pooolt 
ooools 
poe b 
oooolf 
oocols 
oo00is 
00001s 
o0001< 
00no1s 


000015 
o0001* 
000015 
oo0uls 
00002¢ 
oooort 
00002¢ 
000026 


Re aa ik al hl A aid 


i heeeeeeteeeerihatntntiientmemithae hehememnh tenth a anathiieecieninmnnmiaemmamiee ee ee 


TDWORDCLTAD) = 3 
{F (IVALTD.FQ.2) EDWORDCIAD) = 4 


150 CONTENUE 


Cc 


CoccavoeeMASTEREX) ITS THE ARRAY FOR THE REFERENCE DATA WORDS. 


c 


Cacecese CHANGE THE FILE NUMBER(S); THEN RETURN TO CALCULATE THE 


C 
200 


700 


666 
C 


220 


WRITE (6,666) (MASTERIL) sL=1,NPEAK) 
WRITE(6, 700) (ACSI) 2 II= Ly NPEAK) 
FORMAT (LO(5XsF6.2)) 
FORMAT ( LHO, 5X, 513) 


K=0 
Ki = O 
FC = 0 


00 220 IR = 12468 

TF (CPEOWORD(IR)EQ.0)} GNTO 220 

K = K +1 

{( = k 

lf CER.GESMASNL.AND.IR-LEMASN2) 6O TO 225 
th C1C.FQ.1) & = K-ESHIFT 

Y(K) = ArCLD 

ZK) = HUGO) 

6G TO 220 


> KF] = KL + 1 


YOKE = ACKED 
7EKP=AMPLTOCIRD — VCK) 
feo= L 

CON FT NUE 


000029 
000020 
000020 
000020 
000020 
000020 
000021 
000021 
000021 
000021 
000021 
000021 
000021 
000021 
000021 
000022 
000022 
000022 
co0022 
000022 
000022 


000022 
000022 
000023 
000023 
000023 
000023 
000023 
000023 


(+ Oe REEREGRESARERETS KEKE RETRE SHELESSS AS SSESSHS KS OETOESEPSCASS SETHE SSSSE SHE 0002 } 


30uU 
230 


ARTTE (25) TF-Y.Z 
WRITE (23) TF, ¢nWwORD 
IF (51 .GE.NUMBER) GOTO 230 
Jtr = Yt @ | 

NO 309 .1 = IHeNUMBER 
READ (22) TY. FT DOwGRO 
NAETF {255 TV. 1OWORD 
PFAD (24) [ZsAB,HT 
ARETE (25) §27,AB,HT 
CONTINUE 

CONTINUE 

REWIND 22 

REWIND 22 

REWIND 23 

REWIND 24 

PEWENN 25 


NO 310 tL 1=1lyNUMBER 


000023 
000023 
000023 
000024 
000024 
000024 
000024 
000024 
000024 
000034 
000024 
000024 
000024 
000025 
000025 
000025 
000025 
000025 


' 
bus 


Fe ieege te 


& 
5 


SE ES TE OEE A PITT ee ee ee 


READ (23) IV, 1 DWORD 
WRITE (22) TY, FOWORD 
READ (25%) ITT eABsHT 
WRITE (24) TEI,AB,HT 

310 CONTINUE 

Cc 

REWIND 22 
REWIND 23 
REWIND 24 
REWIND 25 


IF (LOOP.GT.MAXX) GO TO 400 


IF = FSPEC(LOF) 
co0P = LooP + 1 
Ga TO 6 

400 CONTINUE 
RETURN 
€ND 
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00007 
00002 
00002 
00002 
00002 
00002 
00002 

00002. 
C0002 

00002: 
00002: 
00002. 
00002. 
00002: 
oncd2: 
00062 

00002 

00002 


OF PETIT TT 


SUBRNUTINE PKTABL CIREG, ITOTAL, IKOUNT) 


90900 


CFC SHeRESREEEEREEEEEE RSE SEEKCREKSESEEEKEEKSERLESSEHEKS £4 KEES ESEKE EEE HEKKER OQDNIAG 


FT TS EL ES OL RE NE PY PG 
aaan 


C 
i 
c 
i ¢. 
my 
‘ 
a 
y 
i 
E 
3 
F 
Pas 
“ is 
110 
i 4 
2 
st 
; C 
3 
C 


230 


ConceccesSEPTEMBER 24, 1979. (UPDATED AT AF/RPL : 


THES SUBROUTINE CROSS-REFERENCES ALL THE PEAK HEIGHTS 
WITH THEIR RESPECTIVE DATA WORD LOCATIONS. 


COMMON /S5PECTR/ TFILEC30) NPT 29(30) 
COMMON /POTNTS/ NWORDS(30) 


DIMENSION AMPL TN(468) » HAVENO( 468) , IOWDORD( 468), AB(S0),HT (50) 


DIMENSION 1030) ,TPRINTC 30) pLEAL 30) 
WRITE PEAK HEIGHT INFO. 


READ C21) TF, AMPLTD,WAVENO 
READ (22) IF »{NWORD 
READ (24) ITF,AB,HT 


WRITE(6,200) LFILE CIF) 
MTA = 0 
NW = NWORDSCIF) 


ON 2 J4=1l_,NwW 

TF (1POWORD(I4).€0.0) GO Tn 2 

MIA=MTA+] 

IF (MIAEQ.26) WRITE{6,200) IFILECIF) 

iF CYOWORD(54).£4.2) GN TO 15 

iF (LOWORD(34).2EQ.3) GA TO 25 

TF CIOWOROCJ4) .FEQ.4) GO TO 110 

WRITE(6,309) "A, JG sWAVENO(J4) FHTCMIA) ,ABUMYAY 
“0 0 4 

WRETECH2 320) MIA, J4,WAVENO(J4) SHT EMEA) »AB( MIA) 
GO TO & 

WRITE(6s 310) MIA, 34; WAVENO(JS4) eHT (MEA) pABEC MIA) 
co YO 4 

WRITE( GH. 359) MIA, S4,WAVEND!J4) ,HTE MIA) ,ABCMIAD 
TPRINT(TF) = J4 

CONT ENGE 


TF CIF.GE-ITOTAL) GOTO 3 
co ro tl 


MXDW= I[PRINT(1) 
10 230 KAN=2,I TOTAL 


{F CIPRINT (KAN) .GE.MXOW) MXDW=IPRINTCKAN) 
CONT TNUE 
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000009 


090000 
000000 
000000 
000000 


CERES KEE TER EEE EE EKER EET EEKESEREKESESSEES SESS ERE HEE FES HBHE KEE KES EH E4EKETENNOOOTD 


ooa0nn, 
onoon} 
On000L 
C0900} 
0000061 
000041 
000001 
000001 
000001 
000091 
ooo0o0na 
000002 
000902 
o00n002 
000002 
090062 
000004 
000cg2 
000057 
000002 
nn00qc? 
000002 
00000: 
00000: 
900002 
000007 
n0000? 
nn0092 
00nnQ 
nn900: 
0g000¢ 
o00nas 
onnnde 
0n07Ne 
onon0¢ 
000004 
000004 
o0nn00e 
oongnse 
900004 
00000: 


‘ 


t 009905 


REWEND 71 000005" 
REWIND 22 HGNAgE 
REWING 24% no990% 

C 900005: 
c on000s: 
C¥eeee = PRINT-ThUT PEAK HEIGHT TABLE. oo0nrs 
¢ on00G®s 
C angnns 
5 1.QUN. = 0 000006 
IKOUNT = 0 ononns 

PD 500 1a=1,20 o00nn«. 
IPRINT(18) = 0 onnnte 

500 CONTINUE 990096: 
ON00Ga! 

on tO Pr=),4xng o000b6, 

DN 35 M=1, 30 onnnn6 
LEAIM) = 0 oneroK 

TUM) = 0.0 coongs: 

35 CANT UNUE ono007 

nN = 0 000007 

T3tGPp = 0 NNrNO7 

45 RFEAN (229 Joy Tl 9 ono0a7 
REAN (24) JdpAAOHT onny07 

IF (TOWORNCIT).NE.0} GO TO 50 000907 

IF (Jg.GFELTINTAL? GN TN 60 nn0097: 

ay 7 45 09096? 

5U N = } 090N67. 

C OnongT 
SPRINVCISs - FORINT( JH #1 ax0005! 

JPL = EPRINT(II) oon002 

T(J) = HTCJPLD onn000P 

1F CIOWORDITE). 5. 27.0R. TDWORDITE)EQ04) LEACII)D = 1 799048 

i oggT Yee 
f. Usd.GE8.FTATAL) 60 TC 60 onc i 

69 60 4 ONU., IR: 

69 FF ON.EQ.3) GO TN 45 090038 

: 90900R 
iy peewee neHRETE “FAK HEIGHTS TO FILC 26 « FOR LATER PRINTING. 00099s 
: 0000G9 
WRITEC 26) Itet unn0es 
TCOUNT = [COUNT + 1 nngnhe 

a o0n0n0y 
Coaccoeeel HFCK TF REGRES WHtl SE CALLED. 090.299. 
c 090009 
IF (IPFG.EQ.0) GOTO 65 990009 
CocceccoeYO PREPARE FOR CORRELATION OF PEAK HEIGHTS AND MECHANICAL 000037: 
Coceccee PROPERTIES 3 SEARCH FOR ONLY THOSE PEAKS THAY ARE VALID IN 0009341 
LwevenceehLL THE FILES. 000034 


c 
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Coesecee Of) NOT YNCLUDE THE WAVNOR PFAK IN SHE REGRESSION ANALYSIS. NAVD 
Ce eawcne sl THE REASON FOR THE NPT29 CHECK BELOW SDETEI+SDIRK) = 0.0). 090934 


TEC EELEG.NPTCICLID GN TO 65 


00 210 KAY = 1L,ITOTAL 900099 

TF (LEACKAYD.EQ.1) GNTN 210 000009 

mrop = t olete kale e 

210 CONTINUE OO0NOLO: 
1© (PSTOP.EQ.1) GOTO 65 NN0914 

C FICE 27 +#¢+4+ VALIN PEAKS IN ALL THE FILES. 000019 
WRITE (27) TI,T oleleres ie! 
CXKOUNT = YKOUNT + 1 000nL4 

&5 CONTING GG 000019 
REWIND 22 erelelon Gen 
REWIND 24 On00010 

tO CONTINUE 9900010 
Ceeee On0GLO: 
REWIND 26 okorered med 

TF CIRFEG.FQ.19) REWIND 27 900011 

C 00001). 
NL - 1 00901? 

N2 = 10 NNNOL?I: 

C 0000) 
50 WRITE (6,410) CIFILECT Ds T=N1LeN2) OONaL): 
MN 240 NI = £,FCOUNT N00011 

PEAN (26) N4,yt NONOLT: 

WRITE €6,400) Nog f TiLK} eK=NL NOD onaNnLT 

TF €N3.FQ-25.9RN3LFQ.50) WRITE(6,410) CEFILECID, IN1,N2) O0n0L? 

e<O CONTENUS 900012 

v 000012 
‘CWIND 26 N6NNL2 

TF (N2.GESTTNTAL) GO TO 70 N0001Z 

TE {IN22EQ2202 “9 10 20 000012 

SS oe 000012 

VR RY 000012 

Gr 1 £-4d ON0901?: 

SO NI = 21 0000}! 

N2 - 365 000n\2 

“fn Th 2730 N00NT3 

‘, 097013 
CER EE MOREE AER EEAE ERE AER EAERERSEESESESE EERESEESOEREEEEORSETEE SH EO4 98S EEO0ND] 3 
C onaaya 
7, WURUTE(6,420) FPTOTAL on00L? 

( 000013 
CORTE ERERES EERELERERRERERERRROERRRES FERRER ERER EERE EERSOSRAEEKE £48 ETE 000013 
C 000014 
{, 000013: 
(&eaee tNRMATS. ONNNLs 
f Q000L4& 


FREE HERETHREKTERELE EERE SRS HAREREKERESEESRASESHE RHE RAKSHA HSRSER ESSE SRE SR RKEONHOILS 


it 


ee ELSIE IT ATT TT TES SPELT LIT PE TOT MPT TE PE PT 


’ 
t 
t 
Hy 
t 
[ 
e 
t 
i 
EB 


000914 


200 FORMAT ELIT S%,"FOURTER TRANSFORM INFRARED SPFCTROSCOPY -- NORMAENODONI4 
1IZFO PEAK HEIGHT INFORMATLON FILE NUMBER®, 13>//5X 000014 
2,*ee — DENNTES A VALID PEAK. "/5X,"CB. — DENATES A SPEAK MEASURED FON00LA 
3ROM A COMAON BASEL ENE 2*//52 000014 
4,"PFAK NUMBER DATA WORD WAVE NUMBER PEAK HEIGHT BAG00014 
SSFLINE AMPLETUDE™/5X_"*—--~-- —----—— eee mare ne ONN014 
() leetelintnelleaiieeetaten Soe ee ~~) 000015 

300 FORMAT ULHOs 8X e129 1 3X0 T3e BXeF lL Oe 4s S5XeF1 2.79 SXeF 12.7} Oo0oLs 

310 FORMATULHOs 4X "*#EF "VEZ gh SX oI Se BX oF LQ Hy 5X yg Fl 2079 5X0 F12.7) 000018 

320 FORMAT LHO, 8X,T2," CRo%eIX9l 39 BX FLO Se SK oF 20 Fe SX yFIL2.7) 000015 

330 FORMAT CLHO 9 4X. BER M1224 CBs IWXe 139 OK FLO. 4s SX oF R20 Ty SX pF i207) OOOKTE 

400 FORMAT {/8X_135 7Xe LOCF 7.39 4X)) o000g1s 


410 FURMAT(LHL, TX,@FT LS *#8 OO See Re eeRe See eeeeeseses PEAK HEIGHT TAAUL090015 
LATION FOR ALL THE FILES 88880 49e SESS RS REE SS CEKEHESEERERETOSERESENDONOL® 


2//BXe "DASA", BX y nn00L® 
31 OCMETLEN, 7X) _ /8X_ "WORD, IXy LOL TL 2e OND / BXy —————— 900915 
Wyn en ee nn nn nn a en 000 
§~--— ~ e+ = /) 000014 


420 FURMATOLHD, 4X» "NOKMAL FND OF PEAK HEIGHT TABLE "/5X,"TOTAL NUMBFROOODLE 
1 OF FILES CISTED = ",1399%2") OON0L4 
RF TURN OOndLE 


END ONOOIS 


t 
t 
& 
E 
{ 
r 
£ 
& 
é 


rupee 


TE SORE MEET E TIE ST TE 


t 
t 
i 
b 


SUBROUTINE PLWAVE ( ETAPE, TF eNPTS, WN, AP) O9)5597) 


CECE CEO CECE CEC ECTECOECEC CECE COCCECCOEELCCECEECEECCECEEECCEEEEECCECOEOOCOEE090 26 29> 


{ 


99005899 
Coeeeeee PINTS AN FNDEVEDUAL FILE SPECTRUM FROM THE FTIS TAPE. 90075904 
( 9000501 % 
Cewe coe ee TUPNKO S OWEN Co. SMITH. -— AUGUST 4, 1978. 0000597% 
u 2NNIVEIA4 
Coanwcen swccccecccccoecce UPNATFED AUGUST 23, 1979. 0N07594¢F 
C 99095956 
COCCE COC CECE CECCE CECE COE CT OOLOCOCOEC CEO ECE CEOEC CEC CECE CECE OCECOCECEEECEEEE 600095946 F 
C 9099597" 
PEMMON /OETME £ FMAX,F MIN: NORM NNKnNSgGAA 
DEMENST OM APLG7O) sWNE47O) » TESTES) 
DNATA TEST/SF IL, NUM® "BER 4/7 10096970 
‘TF - Te 0000692% 
CALE PLDT (CO. Teg *) 
YuIT=e75 CO00N0b9N4E 
YTTTI=A¥TIT-.20 0900Kh05" 
YOU eyTtT7-. 70 OND0bB)45 
YIETGHVYEET3-.70 N00 94076 


CATE SYMAML (.5eYTTTy.L, "AIR FORCE — RPL™,0.,15° 
CATL SYMANE (Loe YTIT2,e1e%FOURTER TRANSFORM INFRARED SPECTRUM%*9. 29007679 ° 


a) ON0 94219 
PAPEL SYHANL CoSgVTETS, ole "TAPE = "yQey 7) NAIVE 22 
CAPE SVM (699, ,999, geleT TAPE 90.94) 00074225 
TAYE SYMBE 1995. 42799 aly ™ Meee?) 00906749 
CAEL SYMBOL (9990. 7990, 01, TEST y 00s 12) DOUIF WE 
CALL SC PIEOL £9990,909. 01 92H= 90092) ONANE PAS 
Crit “UMPER £499, ,969.,e1 FILE »O. 2-1) ON}D%OF F749 
Te ENGWALE OLED CAtt SYMBOL (oS VYTITT4:22le "NORMALIZED SPFCTRAL DATASQONN4IRY 
ie tley 24) OO0NF 29) 
i} (MORMSEDLO) SY TN an 09906299 
YoY MAX 909964149 
MPN OMEN 9900637" 
we Tr att NOCNGHIAAN 

NO0N6340 

SCARING PY oY AXES, 00006359 
ANOAAZEU 

SK AMG) 00096375 

w UME AP OL) 99906781 
1 32 yNPTS 0990632aNn 

1 (APL IPL GESAMAX) = =AMAX=AP(S) 0000640) 
is CAPO SDSL SCAMENY «= AMIN=AP(S) arNNAgI1A 

ro OMT PNUE NANA1 420 

< ATNCP> (AMAX-AMIN) /10.0F0 NI076 +A0 
APINOTS4#]1) = AMIN 0:1)96449 
APPEND ES#2) = ATNCER 0900I5459 

N9N064450 
SCALING MF  X AXES. N000647) 


QANAK4EARA 


, 


b 
& 
% 
F 
4 
7 
£ 
r 
\ 
E 


[ 
E 
; 
r 
E 
E 
F 
r 
E 
F 


SS EE PRR EIT TEE TT PATS STITT tat Aa ak id a ee i eh alla 


MANPYS#Ed) © 42900.0F0 


09074400 
WM(NPTS#2) = -700.9F0 


YRAVIESOD 
TE ENORMLEOL) GO TA 40 NOMIASLA 
CALL AXES €9.,00e¢9NTIRMALTZED AMPLITUDE®, 20, 10.09900 sAP(NPTS#1) -AP (92996524 
INPTS+#2)) 0090965445 
oy 19 50 NANN454C 
oY TALL AXES CO... *AMPL I FUDE™, 9, 10.0990. pAPENPTS#1) APC NPTS42)) 0090465) 
BO CALL AXES €O.sQ. eo" WAVE NUMBER "9-11 52000900 eWNINPTS+#LI YS WNINPTS4#2)) 9000650” 
COiL LUNE UWNOLD ,AP LTD ONPTSo1,0,2) YIINGETA 
CALL PINTE23.0,0.0,-3) 
CR TURN 00004594 
Pt) 


NIIDSE 


SHBROULENE REGRES 00000040 


L 00000076 
t MULTIPLE LINEAR REGRESSION 00000089 
c BY THE STEPWISE METHOD 0000009 
U DEVELOPED FROM THIOKOL PROGRAM NUMBER - £023. 00000109 
c 00000110 
CHORE ERERERE EEE SEER EEEREAAREALESHSERKESRERSEASESSEKSESKE SSS SSSSEEKSEVSC ESTER EOOOOOL ? 0 
C 00000170 
c 00000140 
C N ~ TOTAL NUMBER OF VARIABLES ON INPUT (MAX. = 60). 00000150 
C t - TNATAL NUMBER OF DEPENDENT VARIABLES 000001460 
C M - TOTAL NUMBER OF OBSERVATIONS OOOOOL TL 
‘ Fi — THE F SIGNIFICANCE LEVEL TO ENTER A VARIABLE INTO THE 00000180 
C REGRESSTON 0000019C 
c F2 ~ THE F LEVEt TO REMGVE A VARIABLE FROM THE REGRESSION 00000200 
¢, ite - NUMBER OF TRANSFORMED VARIABLES 00000219 
L J ~ VARTABLE NUMBER BEING TRANSFORMED 00000270 
Cc NTRAN — TRANSFORMATION TYPE CODE 000002 311 
C = Ly LOG OF X -— LOGLOUXEID) 00000244 
C = 25 POWER “TE X - (XCD) +A)@OP 00000? *C 
( = 3, SQUARE RDOT OF X - SQRTEXCIDS 00000263 
is = 4%, NATURAL LOG OF X - LOGEKETDS 00000274 
: = Sy XCI-L) © CXC SP**POWERCIDD 00000286 
{, CONS — TRANSFORMATION CONSTANT IF ANY 00000290 
c POWER - POWER TN WHICH TRANSFORMED VAREABLE MAY BE TAKEN, 00000300 
c 1.0 WHEN NOT INPUT 00000316 
< wo - WEIGH] APPLIED TO THE SPECIFIC ORSERVATION, 00000321 
‘ x - NATA FOR VARTABLES L UP THRU NL 00000336 
> 00000340 
¢, INCORPORATED AL THYOKO MAY 26, 1978. 00000355 
% : 0000036¢ 
i SUARQUTINE “RETIFN GN MAY 224 1979. 0000037u 
¢ 00060380 
RY OWEN ©, 5™- TH, 0000034¢ 

00000406 
Ry Se i RHEE EES SESEEEKSESEKEKES SEE THEKHR EEE EEK EASSRHSS CERES SSOCEEEHOOON0O04IC 
CPMMON Jo RREL /S PRREF(60) »MPFILE(30)0, 0000042" 

‘ NPHYSI (10) sNUMDEP , NUMI ND »NUMOBS, IREG 00000431 
OMUN ¢5tC°TRS LFILEC30) »NPT29€30) 00000444 
YEMENSTON F4YS (100,10) »XHTE30), EVARCS5SODB 0000045¢ 


4 MENSEON FDISTO30)~. SUMX(60),. SUMXX(60260)5 X60), 81603, $8160),0000046° 
1TOINS (60)45D160)5R 460,60) 4 XML 60) » POWER 60) .82 (60) » Y¥ESO) ,CORRCF (10) 0000047" 


>> NTPANL6C) ec FFLAGI609, 1SUB(60) 000004AC 
z 0000049¢ 
Coo wees -f -— OLFSTRERBUTION PERCENTAGE POINTS (ALPHA = .10)6 00000502 
£ 00000535 
Cowes. -- -FRIM THE TEXT » PAGE 482. 0000052d 
C 00000524 
SATA EDEST 2 39086 98.539 505% 0 4e She Se 0s 30 18 e 3059p 3c Shy 303693028, 00000545 
A 3023,3. 1393-6 1493616932079 3005 9 30039 3cOD eg 29D 297, 0000055¢ 
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NN 2 96 9209592 09M 920939 Zu IZ 92 Fb 92 eID e 22899 208992088 / DONDN5EL 


« 
CECH EEE SHEERS EEESREEEREREERESRE RE EHERSELEKA ELSES OS OHS ESSERE ESEHSESEEESEERERERA 


C 
a RFAD MECHANICAL PROPERTIES TAPE 3 GENERATED BY E420 PROGRAM. 
C 
CHEERERRAREKESRE SREREEKE HUES TRS SESLEEKSS EKASRE KE SSHS SF LRASHOHSN NS ORK ESE KEKE EOS 
: c 
c 0000057(¢ 
E REANCIREG) TITLE 
f 4 READ CIREG) NleN29N39N4o (PHYS ENLON5) »N5=1,10) 
I IF (FOFCTREGI) 524 
E N = NUMIND + NUMDEP 00090606 
[ NN = N 000006 !¢ 
: L = NUMDEP 000006?¢ 
E M = NUMOBS 0000063 
t TCURVE = LINEAR GR NON-LINEAR REGRESSION MOD. 0000064( 
CURVE = 0 00000855 
JCURVE = 0 0000066( 
F REWIND IREG 0000067( 
| ltr = 0 0000068 
E {k= te 0000069( 
3 IF (UL.GT.299 STOP 
Fl = FOLSTIUL) 0000071C 
é f2 = FOISTOLLED) 0000072¢ 
{t (1CURVE.FQ.0) GO TO 12 0009073¢ 
“Vos NUML # JCURVE 000007 4¢ 
WRETC Oh, 705) NV 0000075 
omy orm Ly 0000074C 
t WEITF (697108 oo00077¢ 
MRELE (6.720) Neb My FL F2 0000078¢ 
f f 0000079 
t ( INITEACTZE SUMS. 0000080¢ 
t ‘ 000008 1¢ 
7 roy w- ft el 0000082 ¢ 
F N2=H 0000083 
; 1085-6 0000084¢ 
§ HO 20 T=1,N oo000Rs( 
: HAMEED 0000084t 
E HPRAONTT9=06 0000084 < 
, SUAGCTD=0 00000n aC 
t TSUBEL) =0 000008976 
i NN 20 J=LyN 0000090{ 
0 SUMXXLT, 5220.0 9000091 
( 0000092¢ 
NBS=M 000009 3¢ 
¢ tt (TTR) $0,50,40 0000094¢ 
E G 0000095; 
é CALCULATE SUMS, SUMS OF SQUARES AND CROSS PRODUCTS. 0000096¢ 
( RYAl) OBSERVATIONS FOR EACH VARIABLE Xf1). 


090009 7¢ 


ry 
4. 


“ 


wn 


WREIF (6,750) 


() N-N2 
C 
Covee SETUP TNDEPENDFNT VARIABLES. 
C 
te (CFCURVE.EQ.1) GO 10 76 
C 


99 3000 MCN = 1,~NUMIND 
FEAND 27) FAX, XHT 
XOMCOd = KHTCINBS41) 
UVAREMCM) = FAX 

360°) CONTINUE 
REWIND 27 


FIND WHICH RECOKD «. NEEDED FOR PHYSICAL PROPERLY IN®UT. 


t 

C 

C FIEND WHICH RECORD TS NEEDED FOR PHYSICAL PROPERTY INPUT. 
C 


1iR ~ MPF ELECIOBS4 1) 
C 
ComeeeS!t TUP OEPENDENT VARIABLES. 


al 
Y 


00 15 MGM=1,_, NUMDEP 
KIM - NPHYSLONGM) 


X(NUMINO MGM) = PHYSILIB sKIM) 


Ts CONTINUE 

on oro 7e 
T+ PFAD CUS) CXC IAD SI ALIN) 
rae WH "5 0.90 


e 


Canes FPANSFORMATION COSSTIBILEITIES. 


C 
yee 7% Y VELSNn 


1 CHFLA CF)NE.0O} GO TO LT8 


iF CNTRANC LDP 180,180,120 
Po. MINTREN FRAN) 


TO CL 3059140, 1 50016091 709175? eNNTR 


WRETF $6, 780% 


"TOP 

€ 

C ‘rey OF x 

13G KETP=ALOGLOCUXCID) 

Ga Ta 189 

c 

Cc POWFR WF XxX 

C 


149 XCER = OKTL24CONSCL) PO*POWERCT) 
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GOCOVEFA 
00006°9F 
00001600 
00001016 
0000102¢ 
00001030 
0000104¢ 
00001056 
00001060 
00001076 
00001080 
00001020 
o000llac 
oo001lt¢ 
00001126 
000011436 
00001144 
00001120 
00001140 
OOOOLLSO 
0000) 160 
00001176 
000011390 
00001196 
00001280 
000012{2 
0900122¢ 
00001230 
00001240 
00001259 
00001266 
000012706 
o00012K¢0 
00001290 
00001390 
00001316 
00001320 
00001330 
00001 34¢ 
00001356 


00001370 
00001386 
00001350 
00091400 
0001410 
00001420 
00001430 
00001440 
00001440 


Go TN LAO 00001466 


¢ 00001476 

C SQUARE ROOT OF X 00001480 

C 00901490 

150 XCEPSOXCED RES 00001560 

GO 1m 180 00001510 

Cc 000015720 

C NATURAL LOG OF X 00001510 

C 00001546 

160 XULP=ALOGUXEED) 00001550 

GO Ti} 180 00001560 

00CcolS7L 

C XC1-1) MULTIPLIED BY (X€52**#POQWERTID). TRANSFORMATION CARD 00001584 

¢ CONTAINS (A) RESULTANT VARTABLE NUMBER, (8) CODE Sy (C) NUMBER 00001590 

C NE VAREABLE MULTEPLIER (X03)), AND (0) POWER OF X(33- (3 NOT = 1300001400 

¢ 00001610 

170 JJ-CONSOD) 00001620 

XODI=KCI-L IS CXC SID Se POWER CED) 00001639 

GN Ta 180 00001649 

¢ 00001656 

¢ FACTOR TIMES X 00001650 

C 0000167) 

LTS KC EV=XCPP*CONSOD) 00001680 

GN TN 180 00001690 

¢ 00001700G 

¢ VARTABLE NOT CONSIDERED IN THIS PROBLEM. 00001710 

é 00001724 

1TH XtPP=O0.460 00001730 

190 CONTINUF 00001740 

oN 179 T=1—N 00001750 

27 USUBIEDAL 00001740 

' 00001770 

7 CHECK Ti! SFE TF ALL VARIABLES ARE USED IN THIS PROBLEM. (CIFLAG) 00001730 

: 00001799 

eae eae 00001809 

00 18t }=LyN 00001819 

IF CIFLAUUE}2EQ.0) GO TA 181 0000182% 

. {COUNT=ICUUNT#+E 00001839 
Bi a 3oN-T OD001LR40 
OO LAB kI=15K3 00001850 
E KS=E4K2 epades aN 
E X(KS-ECOUNT) = X€K5) 00001870 
E {SUB(K 5- ECOUNTIEKS 00901830 
3 i883 CONTENUE 00001899 
i 1BL CONTINUF 00001990 
b C 00001910 
N=N-FCOUNT 0000192 
4 00001939 

tf tet) 200,190,200 00001940 


190 WT=1.0 00001950 


20u 


CE . 


ars 


(he 
eft 


AAA 


vn a) 


SUM ft GUBSERVATEONS 


bo 210 £-LeN 
SUMR OTD = SUMX CP CORTOXCED) 


SuMS OF SQUARES AND CROSS PRODUCTS STORED BY COLUMNS. 


N0 210 J=1,N 

SUAXXOL, JU-SUMXXOE, S04 (WTR XCEL DEX I)) 
CF CWT~1.0) 27042302220 
OBS=OBS+#WT-1L.0 

TOBS=TAR sel 


“Eft AASERVATLONS GON UNIT Li FOR USE WITH PREDICTIONS. 


to CYOCUR ZECEQLLD GO Lu 235 

WROTE CLPD CXOKD Ss K=L, iD 

WRETE (6,799) TOBS,TFILE( 1 OBS) 
VPREFE (6,800) (XCT3,TSUBCT) »L=L oN) 


te tM -10B5S) 260,269.70 

it LITCURVE.EC.LI GM TO 262 
WRETE (4,310) 

GH TA P6L 

WRETE (6,819) 

CONTINUE 

PEWINI LE 

ict = J 

‘wr O70 Tei 


“Et ANS. 

“4h TV =SUMX CED SNRS 

eAN tt JVEVEATIONS. 

"GbR OMXXOT, DP- COB S®XMCT PE KMCED DD 7 (OBS—1.0) 3490.5 
Cus vel .£4,.1)} 6G TO 263 

‘fF fF. G7.NUMIND) GO TO 265 

ARTTE (6,820) XMCE)D,SOUT? pTSUBCTD,IVARET) 

66 ™ 270 


WPTTE (65819) XMCT) eSDCT) »TSUBCT) 
co Ty 270 


» ACTT= UCT + § 


WRETE €6-¢821) XMCTI5SO00T) pS ISUR( TE)» (PRREFELAD PLAPICT, ICTT) 
icf = {CY + 6 
CONT ENUE 
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00001940 
00001979 
00001980 
00001999 
00002000 
000020'0 
00002021 
00002030 
00002040 
90602059 
000020aC 
00002071 
00002030 
00002090 
0000210" 
00002110 
00002120 
000021 30 
000021%" 
00002150 
00002130 
000021 70 
0000210 
00002120 
0000220uU 
00002210 
00002220 
00002230 
00002240 
90002250 
00002260 
00002279 
00002220 
00002270 
00002300 
00002316 
00002329 
00002330 
00002340 
00002350 
00002360 
00002370 
00002380 
00002390 


00002419 
00002470 
00002430 
0000244 
00007459 


Mt) 300 L-1eN 00002466 
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ROG. P= 1.O 0000747C 

K= 14h 000024RC 

IF (N-K) 300,280,280 0000249 

220 DN 290 J=KeN 00002501 

Pie SD=C(SUMAX (1, JI-COBS*®XMCLDeAM(J))) 2(0BS—1.039/0SDE1 $0055) 0000251. 

290 CONTENUE 0000252( 

300 CONTINUF 000025 

WRITE (6,830) 00002544 

vO 310 T=1,_N 0000255( 

WRITE (6,840) 000025¢.¢ 

410 WRITE (6,859) (ROT,J),ITSUBEL) +1 SUBCI) yJ=T ND 0000257 

WRITE (6-860) 0000258 

a 0000259° 

C CALCULATE FOR NEW OFPENDENT VARIABLE. 00007601 

C 0000261( 

NEA =0 0000262" 

3270 TE (NEND) 330,350, 330 00002631 

330 K-N~L4} 00002641 

TF fYCURVE.EQ.1) GO TO 335 0000265! 

READ (12) CORGT eS) sT2L aN) gp JL SND 0000264 

REWIND 12 000026%: 

GO 1M 337 000076A1 

335 REAT CLSD CORCT SE Sel eND p JL ND 0000265 ( 

REWIND 15 000027 

337 KNOND=K#NE ND 0000271( 

On 340 l=LeN 0000272: 

ROK, EPERCKNEND, ET? 00002731 

340 RCT,KI=RUTyKNEND) 0000274! 

REKPKP=L.O 0000275: 

5N(K)= SOCKNEND) 00002761 

XMik D=XMEKNEND? 000027 TF: 

LO yn 3Iy 0000272: 

‘ 000027 

¢ wt’ ttt POWER HALF GE COEFFICIENT MATRIX. 00002801 

c 000028: 

350 DN 360 F-=1LN 000028): 

Sh a=te] 0000283: 
7 nu, 300 JED 3aN 000284 
Pf re ace © oes | 0000285 
Sou CONTINUE 0000286 
] t 00002R7: 
C WRITE ORIGINAL CORRELATION MATREX ON UNIT L2 « 0000288: 
C TO BF USED WITH SUCCEEDING DEPENDENT VARIABLES IF ANY. 0000289 
E c 0000290: 
c ro RE USED WITH SUCCEEDING DEPENDENT VARIABLES IF ANY. 0000289 
c 0000299 
C 000029! 
; st ({ICURVE.NFL1E GO TO 364 0000297 
WHETC C15) CORED s SD oF=LeN og J=L oN) 000029 3 


SP TUE ET ATTN TEE TOTO RE ETT en MTT TINY 


ATT IE TET TT eT PR PE TE TSE SE IVE Fao 


CRN 8 RT FFE RE OPE FEE ETE BEET TP IE IT BET TS TET TT ET TP | 


AARP RADA 


SNe 


364 


369 
3470 


380 


460 
470 


REWIND 15 
60 (0 359 
WRETE (U2) COKG isd) eT =L Nd eJ=1,N) 
REWIND 12 


REGRESSION ANALYSIS. CALCULATES BEST REGRESSION EQUATION WITH 
STANDARD ERROR OF FACH COEFFICIENT IN THE EQUATION. 


DFGRFES OF FREFOOM 


CONTENUE 
MF =OCS-2L0 
K=N-L4l 


STANDARD FRROR OF OFP.L VAR. 


NVIN=U0 

NVOUT=0 

NTRY=0 

SF=SN(K) #1 (088-1. 0) *R(K Kd /DE 48.5 
NTRY=NYIRVELE 

K=K-1 

VMIN= 1.0635 

VMAX=0.0 

NMIN=0 

NMAX=0 

on 470 f=1,K 

TF (7° CL, 19) 390¢470,390 

1F (ROE ,F3-0.901) 470,470,400 
VPNCHR (ELE, KHL) FRUK et EP/RET ST) 
if {VRNCE 440,470,410 

fF (VRNC-VMAX) 430,430,420 
VMAX=VRANC 

MAX] 

B17, 

SUE: O00 

Gry Tir Gyo 


ALCULAT! R-TGRESSTON COEFFICIENTS AND STANDARD ERRORS. 


oo ACAP-KOLS KALI CSOCK ELI /SSOCEDD 
SP CLI=HCSESSOCLII® CREE 1) 71 O0BS—-1L.0) )@#.5 


TF (VMIN)D 460, 450, 450 


7 VMIN=VRAC 


NMUN=E 

GN TO 470 

li (VRNC-VMIN} 470,470,450 
CONTINUE 

SU4BT=000 

) 480 T=19K 
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00002040 
00002°59 
00n0. 46 
o000z¢e7¢ 
000029480 
00002990 
000030010 
00003019 
00003076 
00003070 
00003040 
00003056 
00003040 
00003076 
00003080 
00003090 
00003100 
000031!6 
90003120 
00003130 
00003140 
00003140 
00003 1uot 
000031 TQ 
00003180 
00003190 
00003200 
00003219 
00003220 
00003730 
00003240 
00003250 
00003250 
000032 70 
60003230 
00003730 
0000336U0 
00003310 
0000337) 
00003330 
00003340 
00003350 
00003350 
00003370 
00003320 
0000339) 
0000340 
00003410 
0000342) 
00003430 


BAe 


’ 


~ 
. 


aanan 


Se was 


SIE yy 


480 


4Al 
482 
490 


500 
510 


5290 


530 


540 


550 


54éD 


570 


SUMBE=SUMBEtBOE)® X01) 
REGRESSION (PURE) CONSTANT 
BO=KXMCK +1) -SUMBI 


OUTPUT REGRESSFION COEFFICIENTS 
STO. ERROR OF COEFFICIENTS 
STN. ERROR OF DEPENDENT VAR. 


TF CNFRY.NE.1L) GO TO 481 

WRITE (6871) NTRY,KSC 

GN TO 482 

WRITE (6,870) NYRY 

IF (NVOUT) 490,500,490 

WRITE (6,880) ISUBCNVOUT) 

WRITE (6,890) FoUT 

cn TO $20 

JF (NVIND 520,520,510 

WRITE €6,900) YSUBCNVIN) 

WRITF (6,890) FIN 

WRITE (6,910) BO 

WRITE £62920) (BCT) TSUBCT) 1 =2_K) 
WRITE (60930) (SOC TI), T=1,K) 

WREYE (62940) SE 

(F (VMIN) 530,550,550 

[tf (VMEN*DF/R(K41, K+1)4+F2) 550, 550,540 


SFLECT PIVOT ELEMENT FOR MATRIX ENVERSION. (REMOVE VARIABLE) 


*¥ 0: NMIN 

NVOUT=NMIN 

NVIN=0O 

TOUT= (VMINFDF/R(h41,KO1) 9) *€-1.0) 

Ne -OF #1.0 

GO TM 57) 

TF (VMAX4OF/(R(K4+1,Kt1)--VMAXI-FL) 660,660,560 


(Ain) VARTABLE) 


KO=NMAX 

NVIN=NMAX 

NVvOUT=0 
FIN=VMAX®DF/(ROK¢1L,K*25—VMAX) 
NF -NF-1.0 

K=Kel 


CALCULATE NEXT MATREX. 


TON FOR NEW NEPENDENT VARIABLE. 
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INVERT ONE VECTOR, STORE ON ORIGINAL MAT-— 
RiX. MUST READ BACK IN ORIGINAL MATRIX WHEN SOLVING NEW REGRESS- 


00003440 
0000345¢ 
0000340 
00003470 
00003486 
0000349C 
0000300 
00003516 
00003520 
00003530 
06003540 
00003550 
00003560 
00003579 
006035aU 
0000350 
00003600 
00003610 
00003620 
00003630 
00003646 
00603650 
00003660 
00003670 
00003689 
00003690 
00003790 
00003710 
0000372C 
00003730 
00003740 
0000375¢c 
00003760 
n00C3T7O 
00003780 
00003790 
00003800 
000038; 6 
00003826 
00003836 
00003840 
00003850 
00003860 
00003870 
00003880 
00003890 
00003900 
00003910 
00003920 
00003939 


Aaa 


AMO 


AAO 


Rap eres Seats 


¢ 


TAS 


5RO 


590 
600 
610 


620 
630 


640 
650 


669 


66} 


00 o10 t=1,K 
TF (I-KP) 5905610, 580 
00 600 J=i,Kk 
IF (J-KP) 590,600,590 


E NOY KP, J NOT =KP 

ROC JP=ACREKP KPDERI Ty SI-RET eKPP OR (KP, 5D) /REKP KP) 
CONT TINUE 

CONT TE NUE 

N90 670 T=1.K 

1F CT-KP) 620,630; 620 

1 NOT =KP, y =KP 

RE tyKP)=(-ROET KP) /REKP KP) 
CONT ENUF 

Of 650 J=1,% 

TF (J-KP) 640,650,649 

1 KP, J NOT =KP 


RKP, SI=ROKP,ID/REKFAKPD 
CONTINUE 


' =KP, =KP 


RIKPeKP)- LO/RIKP KP) 
50 TO 380 


FESOILVE FOR NER DEPENDENT VAR. 
HENI=NE ND 41 
“U2 T7Lt CORRELATION COEFFICIENT 
wen te} 
PRE TofE. Od [NF &5Fee207((0BS—-1 OF FSD(K)SS2 ) ae.5 
ei.Tf 16,950) RMULT 
tf (ICURVE.EQ.1) GN TO 661 
OR ENT LAS? SET OF COEFFICIENTS FOR THIS RMULT. 
wRETE (23) (BCTAD),1AD=1 ,NUMI ND) 
WRITE(6,920) (BEIAD), ISUB{( IAD), EAD=1 2NUMIND) 
CORRCF(NEND) = RMULT 


GQ TO 666 
TH CRMULTOLT.CORRCFEIICURVE)) GO TO 665 
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00003740 
00003950 
00003960 
00003970 
00003920 
00003990 
00004006 
00004010 
00004026 
00004030 
00004040 
00004050 
00004066 
00004070 
00004080 
00004090 
00004100 
00004119 
00004126 
000041 *0 
00004140 
006041506 
000041460 
00004170 
06004180 
00004190 
00004200 
00004219 
00004770 
000042~0 
00004240 
00054250 
00004260 
00004770 
00004286 
00004290 
00004300 
00004310 
00004320 
00004330 
00004340 
00004350 
00004369 
000043706 
00004360 
00004396 
00004406 
00004416 
00004420 
00004436 


665 
Cc 
C 
c 
666 
C 
685 
686 
P £89 
| 
E t 
3 Ct 
c 
a C 
E t 
E 
: 110n 
Ll 
C 
C 
1209 


Ee 


WRETE(6,72000) NV 
GO T0 666 
WRITE(6,1010) 


PREDICTION OF DEPENDENT VARIABLE. 


WRITE (6,960) KSC 
YPRED=0.0 

SUMD = 0. 

SUMD2 = 0. 
KPRED=N-L *#NEND 
KYY=AN-L 


REWING 14 

NO OBO T=l eM 

{F (CICURVE-NE.L) GO TO 685 
READ €14) (X(K} gK=1L NI 

GO TU 686 

READ (11) (XK), K=1,N) 
CONTINUE 

O00 470 JS-leKVY 
YVPRED=YPRFE DEC RC IVEXC ID) 
CONTINUF 

YPRED=YPREN+B0 
AXX=X(KPRED) 

NEL TA=AXX-YPRED 

SUM) = SUMN + AXX 

SUMD 2 = SUMD2Z 4+€CAXX)*#2 
WRETF (6¢970} I,X{KPRED),YPRED,DELTA 


YPREN=9.0 
REWIND 11) 
RFWIND 14 
NeyY = VARTANCE (FROM PAGE 210 IN TEXT). 
SNUMR = UMD 2 — ( (SUMD) *##2/) 


NEVV .- SNUMR 7 (M1) 
S1i1GMS = STANDARD DEVIATION. 
“TOxMa= SQRTINE VV) 
KBARR = [SUAD/M) 
CVV = (CSIGMA/XRARR) 
ARETE (£6,980) XBARR,SIGMAsCV¥ 
KSC = KSC #+ 1 
{f ({t -NEND) 1200,1200, £100 
WRIFE (0,860) 
on TA 320 » 


START CHECK FOR NON-LINEAR REGRESSION ATTEMPTS. 


TF {kCURVF.EQ21) GO TO 1210 


00004440 
00004450 
06904469 
0000447C 
000044560 
90004470 
0000455ft 
00004510 
00094526 
00004530 
00004540 
00004550 
00004566 
00004570 
000045H0 
00004596 
900460:0 
10004610 
00004626 
000046 34 
00004640 
00004650 
00004660 
00004670 
00004689 
00004690 
00004760 
n0004710 
00004726 
00904770 
000047T4C 
0000475G 
00904760 
00004770 
00004786 
00004790 
60604800 
00004810 
00004825 
00004830 
0000484C 
00004859 
00004860 
00004876 
00004880 
G0004890 
90004990 
00004910 
00004926 
90094930 


Rerewee ans z 
OS ET aE re a x rue ry en 


EET 


pe 


elit All at bo ee Td 


aS hte ras 
TEENS Oe RRM CMMI HONS 


1210 


1250 


9999 


AA 


1310 


70% eORMAT (1H1,44X, 


" TY FORMAT CLH1944X,"E690 ~- FTES REGRESSION ANALYSIS"//1H ) 


720 FORMAT 19H NO. 
1! ES yh3//22H NOW 


TAY FORMAT 


ENTERING NONLINEAR REGRESSION MODE. 


REWIND 123 
ICURVE = 1 
MIAME = L 


NUM] = NUMIND 
JCURVE = JCURVE # 1 
IF (JCURVE.GT.MIAMI) GO TO 1500 
IF (CORRCF (JCURVE) .LT..750E0) GO TO 1350 
READ (13) (81053)5J33=1,NUM1) 
READ CLLD CVYCEED S TI=1,_NN) 
TF CEOFL ILD) 1400, 9999 


iA = 1 

DO 1310 J4 
IF (B81(J54) 
X€TA) = YY 
X(HTA41L) = 

IA = JA + 
CONT TINUE 


=, NUML 
«FQ. 0.0E0) 
(54) 

XC LA} #*2 

2 


GO TO 1310 


X€TA) = YYONUM) -JCURVE) 
WRITE (14) (XCKA) eKA=1,I0A) 


GO TO 1250 
READ (13) 
GO TO 1210 


CONTINUE 

REWIND LL 
REWIND 14 
N TA 

L L 

50 TO 10 


tou 


FORMATS. 


(ICURVE=1}. 


V/SRX,"UCPENDENT VAREABLE = VAR(®,12,°)°//1H } 


"SLE 4 F1003,34H 


OF VARTABLES 

OF OBSERVATIONS 
F LEVEL TO REMOVE VARIABLE 

CT2eT2eFl004,F1 0.4%, SX,EL) 

14) FORMAT (9H VARTABLE,I¢,20H TRANSFORMED. 


745 FORMAT (9H VARTABLE»1I4,20H TRANSFORMED. 
X> THIS VARTABLE 1S NOT CONSIDERED PART OF PROBLEM ",19/1H } 


2 13,34H 
015//28H F LEVEL TO ENTER VAR 1A00005360 


eF10.3/1H ) 


00004940 
00004956 
000049640 
000049707 
00004980 
00004990 
00005000 
00005010 
00005020 


00005056 
00005066 
00005070 
00695080 
90005096 
00005100 
00005110 
00005120 
00005130 
00005140 
00005150 
00005160 
00005170 
000051380 
00005190 
00005200 
00005216 
00005220 
00005230 
00005240 
(0005250 


oeete FTIS — NONLINEAR REGRESSION MODE FOR **#8*%"90005260 


00005270 
00005780 


NG. OF DEPENDENT VARI AB00005290 


TYPE, 14/1H ) 
TYPE oe 1%s3X-" 


00005310 
00005320 
00005330 


® IFLAG=100005340 


00005259 


Te O FORMAT (LH ¢24HTRANSFORMED DATA VALUES .//1X,"OBSERWATION / FTIS F100005349 


T3830 FORMAT 
790 EQGQRMAT (1H 
430 FORMAY 
Lo y1lH »F10.05¢5H VARE,12, 3H) 


tLEONOL"7/) 


{7TH 


POH VAR (212534) ) 


eF10.5,5H VAR( 212 43H) 
el »FL10.5s5H VAR( 12, 3H) 
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{44HCOMPUTED GO TO INDEX OUTSIDE ALLOWABLE RANGE) 
oT4e"7% 512) 


00005370 
00005380 
00005390 


vlH «F10.5, 00005410 


00005420 


810 FORMAT (LHL/1X,"* MEAN STD. DEV. VARIABLE 000054 


1 DATA WORD / PHYS. PROP.*/) 000054 
B15 FORMAT CLHL/1X," MEAN STO.~. DEV. VARIABLE 000054 

1*/) 000054 
BLO FORMAY (OXF 10. 4,4XE1L 2.4% 7X16) 000054 
820 FORMAT (9XF10.4,4XE1L 2.4, X16, 16X,13) 000054 
B21 FORMAT (IOXF 10.4 e4XEL 2045 TX 1691 6X9 6A4) 000054 
830 FORMAT (1H /47H SIMPLE CORRELATIGN COEFFICIENTS. (ROW BY COL.3/1H 000055 

1) 060055 
B40 FORMAT (1H ) 000055 
BSO FORMAT (1H F72.4s4H RET 2, tHe el 2e6H) FT e4eSH Rig 12,1 Hes 125£H000055 

i) oF744e4H Ry T2,1He e125 6H) oF Te SeGH Ro eI 2ehHee ¥206H) 000055 

2 oF %e4e4H ROG 12s LHe eT2e6H) ) 000055 
B60 FORMAT (€1H1) 000055 
R70 FORMAT (1H /14H TRIAL NUMBER ,I5/1H ) 000055 
B87L FORMAT (1H /14H TRIAL NUMBER 915e3Xe"FOR VARIABLE €°,12,°%3"/1H 9 000055 
880 FORMAT (23H VARTABLE GOING OUT = 4f4/1H ) 000055 
890 FORMAT (20H F tFVEL +F12.4/1H } 900056 
900 FORMAT (23H VARIABLE GOING IN = ,I4/1H ) N00056 
9LO FORMAT (20H PURE CONST. BOO) = ,€12-4/1H } 000056 


920 FORMAT CLX"CUEFFICTENTS*//C1XEL2.4%5° BE, 125%) #,E12.4%," AC, 12,000056 
bd F124 BIG e F2_%) MeEL2e4e™ BEM T2e™) @2E12.4," BI",12, 000056 


2m" )) 000056 
930 FORMAT (1H /31LH STANDARD ERROR OF COEFFICIENTS // 000054 
LELXEL2 049 TXELZ 04e TXELZ 04s IXELZ 045 TXELZ 049) O000SE 
940 FORMAT (1H /28H STANDARD ERROR OF ESTIMATE //TH 2E12.47/1H 1000056 
950 FORMAT (1H //33H MULTIPLE CORRELATION COEFFICIENT//7H ¢F10.5900056 
4) 000057 
960 FORMAT(LHL, LOX,"ACTUAL VS. PREDICTED RESULTS FOR VARIABLE (%,12,%)000057 
I*//1X,*ORSERVATION™ 5 BX, 000057 
PM ACTUAL PREDICTED DEVIATION *//) 000051 
970 FORMAT (175F20.495F16.4/1H } 000057 
DAD FORMATI/IXs*MFAN = % pF 1206, 4Xe"*STDe DEV. = "eF 2065 4X_ "COEFFICIENT O0005S7 
1 OF VARTATION = ",F12.6)} 000651 


1000 FORMAT( LX, **e* CONCLUSION NON-LINEAR REGRESSION MODEL PROVIDES 900005) 
LTHE BEST FIT FOR VAREFABLE (7%, 12,%) .*/1X, MHC RSSESHSREEHESHESESE4EE#D00057 
PHESKEEEKERE SEEKSRAE FT HEKEEKKERERS KSRSERERSSSES SE OS VOSSSHSSSS EL SSSS EE SH00005) 
jae) OOO0S5E 

10lU FORMAT( EX, *#@#® CONCLUSION LINEAR REGRESSION MODEL PROVIDES THE 000055 
LT FIT FOR THIS VARIABLE, S&@e® J 1X, FRASRAEKSS SHKSHSSSSSSSSSSSHS YS 4 000055 


PRERESHSPERKEKEVERESHESSSEASRSESHKREKSERSES EAS ASS SS SCR SSHRSe ) O000SE 
1500 CONTINUE 000058 
RETURN OOOO0SE 


FND 00005F 


PRAGUAM TALTAGTMEPYNPUT BUTT TAPS O: ENCUT, PAP FEE NUT PUT, TAPES, TAPE9) 
COVERS REED KEES OREO ED REE ETRE E REEL SERRE BER ESS EERE RECHERERS TERE EREEAE DEE 


C 

C E41U 2 FUR dee TRANSEOAM TNR AREY SPEC TRAOSUCIY. 

C 

C OHYTTCAL PROPFRTTIES MASTFR TAPE GENERATOR, 
c 


CCH HLSRHE KL OLY EERE T CPOE EC RRERAEEE AES ESERARE OE ERES SRN KEKE S EEETE TE 


THENK A COOPORATION 7 HUNTSVILLE , ALARAMA 352072 


PRINCTVAL TAVE STIGATOR We We SCHWARZ 
TELEPHONE (295) - RRZ2 — FINK 


SCTEMTTETC Sr OAGRAMMER ON. C. SMITH 
TELFRHIONE (205) — RAP — R215 


RMARHELRTHR ERASER OTL LET EKER CE RNG TAKEO SHRTS SEKI S CHREREEERETRE TER EEE EES 
SFPTE MATER 4 47 79, 
Pear aM TV - R LEVEL OL ' aNGJAGE 
COP BAN == TAFRPL) 
PREOALC UV TM UL PILLYES T CF CONTHACT FOGAIL = 78 - © = 0027 


PELE LISP LARS EER EL GER RRA HERES SAREE SED RSL EH AHR KEKESE SEWER EE 


BorsaAaaAaanrerraanananwarnanan 


EMENE TAN MECLG LOD TITET CUD)» TAGEWK IL OND, TTEMPACLOOD,ITESTT COLON) 
SPMENSTAN PSEC ETON) 

COMMON NL YS? 

NATA CREATES 8OCO ATE FMS SUPHATE/PUPOATE "/ 


C 
¢ 

MAS Ter = 8 

SFo= § 

nN? = & 

Peo 

REMY EMEC D Cotaree 
C 


TPE Te UNTYE SCORES YATE) GO TA 480 
fF ERUNTYS MFCR CATED GA TP 20 
GM Te 3 

> CUNT INGE 
WRETEIN2, 2 70) saInpTyo 
SUNP 71 

4 CDA TINGE 
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PAD O%1,479) TITLES 
WRETE (N2,S270) VITLE 
wPTTPON2, 5469) 
WRITE €08:7,530) 
WRITE (N2_,%%5) 
WRITE UN2, 540) 
WRITE (N2,5590) 
WRETF (N29545) 
WRITE (N2,559) 


TCOUNT=121 
LO VEAND ENT SIND FLAG ge KL K2 pK? 

TE (EOE CML 200,096 

999 TCOUNT: PC nUst +? 
NREGSNREC+T 
TAGEWKENRES) = KI 
ITEMPACNREO) © Ke 
CTESTTON@QES) = K3 
VEAN €N1,510) CPHEINREC 4), ‘=b,19) 
PRETF UN? 200) SoEC,PAGEWKUNREC ), TTEMPAC NRE) QTTESTT(NREIT), CPHINRE 

—Cy SJ), I=T,1O) 

TREC ENS FC DSNREC 
PE QCIC CUNT SE CAS) Gst TN 49 
Gn Te 10 

40 CONTUNGE 
YRETECNI, 579) 
WRUTE (N22, 849) 
Ww? ETE (ND ,S8AN) 
WE ETE CNI_635) 
PETE CME) 
AQIT! CMP, ey 
MP ETE (2,545) 
MOETE (I, SEN) 
CCM = 7 
ero or ya 

: 

650 CONTINUES 

CATE OF NAME 


heoT 260 


AYO CUOMT EME 
WOTTE PMASTOAD PETES 
Ww) 246 N= } ‘ BST r a 
SWPETO OMASEE OD PIECONDs FAGEMK OND, ITEMOAIND,S TTESTYIN) gC PH(N, JUD, 
LIJ=L, 19) 
247) CONTINUE 
C Pa 
CORRS ERE CEE ERRERECEE REE EEEEREREERERRESRERERERRERTE SE SHERMER AEEE KS SOE YE OE 


C 
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479 LOMMATELIE SER, "C0 AR," 1S AN ENVALTN sunf. NN CHANGFS APPLIFN T 
JP.) DATA-ShE. © XOTCUTION TERMINAT TON eS ° oo 7, 1th) 

GWO COMMAT CLVA4) 

4QD VERMAT CAR) 

BOO ( IRMATCAR, Past, 2010) 

S10 FORMAT (SCF 10070) 

70 FORMAT CIMDSEYS@F WOLTER TRANSFORM ETNERAREN SCFCTRASCNPY ~ FSi ote 
-YSICAL PROUPFRITECS MASTER TAPE GENERATOR" 
WPI UMEAFCPL ZOO -— THINK ZICINTSYVTILDUE DT CZUXPPRNPERTIFS Fo? STLIN © 
URNOELLAMNTM sD SK eT DAS? : 

SAO FORMAT (2X,"REF AGF AGF TEST MNMULUS STRATN AT STRAIN AT 


~4MAX PMs Sie aim STRAING) 

£45 PORMAYT CEN, ONE, TEMP TEMP BREAK A KTM 
— STRESS ENFPGY FANURANCE®) 

S40 EONPMAT (450, 9S7T~ FSS6%, 15X%_ MNENSIT TY) 

545 EOPMAT (5%, 7CWK) CF) (=) CPST) qt) tt) tac 
-1) sey) cc2@) 

SSN FIYRMAT CX, Man - ne - eee tenn - tee nee en ee. ee ee eee mee - + ee ene 


ae ae | 
PVD FIOMAT OZER Ee Whe LHe Le TH, ORT te DK OL OLE DF. 3, 779d 


Rep COOMAT (1 0, Sm ot + - tere cae BETTE ARASee Ses SVE gas :TERABRAART St ge Hes -2> 


t- 


rT LT SV PMS Le CMS SE SRERS SRABRE MoS UMRSEC LK*TK Es TE FSRSVE CFRESARURTBKit cS 


; Ys seere22@) 
SPI EMRMATV GUN) 
c 
CHLEI*SNHHESRSESRAESHHS SHAS HHHRSERSSEHRAESEHRHS EHSEKAHAHSSE CHS SSHSS ARS SRS FSHUSSRSE FE 
¢ 
2990) CONTINUE 
OP MPN MARTE? 
NTA 
rua 
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SUBPCUT INE RENA Sr 
LECHER AHS AEE SE SERAHERERESSEHSSSHERARTSHSESSLHSESE TEKS SESHESHKAEATE SEH ERE 
C 
Ceccacecaen ot PTE Mb YD 9 979, 


Aan 


CHANGES THE OATA SFT FSTAAL ISHED IN ©6410 MAES, 


sy 


CESHHHHSHOHS SEC SHSTSERSOKASLSHSHH ESS SHCHCSHSSSS SS SES KS SSHS SS SSRSSASHRKSKRLSSESHSEE 
COMMON AIT, 8? 
NEMENS EER JICEWK (LOND, STE MPAC LOD) AUYFSTT CLO) oNREC (LOO), 
TPE Ted, 190, SAVOA TI LON, 19), TAGENK( LOD), FTEMPAL 100) TTESTT (109) 
VEMMSEON KOUNT CIO) VILOVETUTLE CIO) 
VATA CRDEST APF NIL ESTES CHANGEJ*®CHANGE fF DELETE/EDELETE 7, 
LADOZANN #7.7 NONVL J PENDNOL STe7 
MASTOR <8 


WHITE ES 7, I) 
LIntnrsco . 
t-.] > 
PTAYD CCASTER)D TITLE 
1? Pe ADUMASTE RZ) NEO CT DV  SAGTURIT ITE MPACED OSTESTTURD oC PHIC? JU) psa! 
—~21) 
TE CECE OCMASTERDD Leg 99 
YR TEENO PE CED GTC INTATD) CENT OTANRECEED 
1-4} 
ah Te Yo 
Ur CONT Peg 
SOP WEA MART Ce 
F PHA SHH TARAS RS PE GHCSKC SS LHKERSEHKSS CSR SL SHS SS LEKSSAETS SHSASSRH SSS SN SEES YP OE 


C CHEFS Tee TMMYT STONENCE COR ALL THE ACTION CARNS, 
¢ 
CPTHe HSC VEAR FHSS SK KHS AKG SHHSRSSHHKRSSSRSESS ESSE FHRESKHS SEE SHGRASSTSKSESS RAVER EERE 
tT YRS =z 4 
ISHUFT = 4 
LAST = 4 
ar hy 
+i ru 
“tye Ff = 6 


S REANINT, FQDACT Eat NE RD KI ,K2 

T© CENEENTDD Worry 

PY | erwo = CNN e& | “ 
TR GRE UEP ME SCUANGE CANOE OWPLEQLLD PFIRST = 1 
Th LACTISNCE I CANGE) EA YS 70 
TE CAPTTIAN.F S.A) GO TON «EF 
TE UAC TEOMFQSVLETED LAST = 1 
fr Ey a4 


PO TEVER TE STFOL TANI. LONO GT .L) Cm Tn 65 
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fe UMTOOEF ONT C1) GA TO 26 
IF (UAST.YFLY) GF OTA 25 
WRITE (6,710) LONP,ZACTION 


STGP 460 

55 WRITE (6,719) LOOP,ACTION 
STOP 62 

AO MIPMF = | 
Li=tiel 


KOUNT(LT) = N@ 
TF (LAST.E2.0) GF TH 256 
WAETE (6,710) ENNOZACTTIONM 


STAP 64 
WPYTE Ft vEMP, FSK FILE 9. 


Aon 


25 APETTE (0) ACTLON,NR KL KI K 
Tr CACTUS AAT. ANDSACTIV NE CHANGE) 60 TH § 
VEAVIND S 275) CYOKD 9 K=1,10) 
WRIT? (7) (VIK),K=1-10) 
Gn TVs 


30 CunTrruies 
PEREND * 


ed 


LS PrAN(D) “SCTEON UR KL K7 AX? 
(FO CEOTES)) 65.,4666 
7 566 TELACTION. FOE YONOL) 60 TO 690 
TF CACTIOMLFOLE OE ST) GN TO 500 
TPONK UF. TD) Gn ™T™9 450 
TEENS CGT.) Gr T9 4509 
TECACTIONSEOQ CHA CE GA TO 1909 
FR CACTESNLEQLIFLETF) GD Th 200 
TEC ACTION. T 9. AD} 6n Ta 390 
Gy TA 400 
FARR NA RHETT ERE RARE REE FE ERRERKEKSESREEERS KERR SEER E ERAKAKEAERECAAREAE KEARSE 


7 ORUTTON = CMANGE ROUTINE 
CERREREREE ELE CEKEEKL TRE RRYEERIREA TERE BEREREESESES & CE SESREREEEEECERESERD 
' 200) CONTENGS ; 
SAGE WK INP S =#) 
ATE 4ONQENG 2K? 
JTESTTONS P=K2 


V1 MnQeo Syms eee Em ENORER cen LL 


} 

E 

| c 

PFADTAD CPMLIND, £),F=1,19) 

: C 

: . WOETE(N7,7 35) Me 

: WRITE Ce L7IDNK. SAGEMKENR Dg STEMPAC NR » JTESTT INR DS SPHLINR, JUD 
; —e SI=be19) 

: TEGMRGTLINTCT) CYNTOT=ENR 


" CA TR YS 
{REEL REREH TERE EREE ERE EEAERE ESS EKERE SES EEAKRERREE ERE EORHREE SEE SEEKS ES F & 


APNE RED Oar orcas nnanan.teupemerenwmenee Concen o - 


C ACTION = DELETE POUTINE 
| CORSE ED C42 4OEEE SEER ERESSETED SES OE SOS FEES SESEEE £OESESORESOREREEHEEERE 
200 CONTENU 
00 205 tX=1,01 
TE CM2 STK UINTELKDY ESHEIFT<ISHIFT 41 
205 CHNTINGE 
MI = NRISHIET 
NO 230 J=MI,LINTOT 
= J+1 
TFQM.GE.100) STNP 7? 
NR FCC ID=ENREC( 4) 
JAGE WK ( JD= JAGE WK 4) 
JTEMPALSINITFEMPACM) 
JTESTTCSV=5T*STTOM) 


AF 729 L=b,19 
PHIC Ss FP PSPHIEMSE) 
220) FANTINIE 
230 CONTINUES 
CUINTOTHALINTV-1 
WRYITEIN?, 245) Re 
ao TH, 15 
COEREEERESEEEEEKEE ERED SERAEREKE SE CERKET EF F5FERHERE SCHRESRSEOREKEEERAKEEBRE EE 
C ACTION = 4D POUTINTS 
CHEHSEAELEEHESSETCAESE SEES FEES AAESERESE RAE FESS EREERS O8SEEOESEERETEEREEESREE 
3090) COMTIMUL 
Cc 
TE CRP CT.LINTOT) GO TO 340 
AA B2Q J=LslLiInray 
PFEMP FOLIPEC( S35. SP TO 41 


Gey 25y 


490 CORD TSS 
PSL HING Se) 
i" Ar) Kao Bet 
MW; 3L4 t=1,10 
SAYDAT(X PDH CK, 1) 
ee WON TENE 
USGEWK (KD = SAGE WK OK) 
UTP MPALK PZ JITEMPAL Ky 
GTESTTV CW PSSTEST TOK) 
320 CANTINGE 


Mme’ 315 K=5ek 

a oe | 

ay Ws tey,ie 

PHLOEM, EDSSAVIATOK, ED 
ALP CONTINUE 

JAGFRUK EM) = TAGE BKK) 
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eerie 


{TEMPA(M)=TTEMPA(K) 
STESTTU*PHTTESTTI KD) 
CONTINUE 
CONT ENUF 


CINTOTEL 
AN TN 3590 


CONTINUE 


LINTOT=LINTOVeL 
CONTINUE 


JAGEWK (NG 5 =K 1 
JTEMPACNS D=K2 
ITESTTON®)HX3 


READ(O) CPHTINR, 1), f=1,00) 
{F INP .GTALINTIT) GOH TO 3280 
ARTTFUN[], 740) NR 
qn 1s 

380 BRETEI(ND,5ES) NR 
6H Te 45 


C 
COARSE ARSE TERETE SEER EEKEESERKAHESLEERKERE REN SERS SE SHSESESESSEESECRKESEKSSEF F 


© ERROR MESSALSS GP TNVALED UPDATE ACTION AND BAD CINE NUMBER 
: CORES ERY TERETE ER EF SE REEESEESEREE SER EERE ERS RECESSES BE SESERE SE EERED EE EES ET 


GAD CON TINUF 
WUITEIMN> SIO) ACTION 
WRETEON2, G70) 
an TA Sag 
a) CONTE NUE 
WRITE (N22, 469) “8 


stor 7% 
CRA RE EYES TERE RHEE ERE SEEREERESEEEESEEERESHESEERES SERRPREEEEEREEEESEREE 


C CCT EMMSE’OPST TLINES RENUMAERED. DATA STORED AND PRINTED) 
CERRA ELSE EES EREREEERESEEERES FHREEERERRE ERRERESEREFESEROEE EKELE ETERESHS 


Rae WOT CHASTE) TITLE 


EF re ee Ee ee eT 


: “hLTigih?,39)) THTILE 
“PETE (N2 9579) 
WITTE (N29 S20) 
WRITEINZ 2 540) 
WRITE (N29 £59) 
WP ETEON2 (549) 
WHITE ON2, 565) 
WRITE (N2, 560) 


steeergttes 6. 


“ 


1CNUNT=AL1 
AN 58D J=L,LINTOT 


HRD, MMA RRARRARRIOIIN <o e sagtsemeneager sean 


i 


IC QUNT= TCOUNT4 2 

NPECQSV=S 

WRETEC(N7,519) NPEC (5) 9 JAGEWKC(S) gp JTEMPAC SDC UTESTTEO UD: (PHL Jed De I= 
11,10) 

WRITE (MASTEQ) MOEC EDD» BAGEWK(S be UTEMPACS PQ UTESTT OSD (PHC ST) I= 
L1,10) 


IF CICOUNT.GF.AL GN TO 475 
Ga tT 530 
§75 CONTINUE 
WRETE(N2, 501) 
WRETECN?, 57°) 
WHITE IN? 4570) 
WPITFIN2, 540) 
WRITE {N2,559) 
WRITEIN2, 560) 
WRITE(N2, 565) 
WRITE (N29580) 
ICQUNT=7 
$80) CONTIMNUF 
GN Tn 450 
LEREELEREREREHRERRESEHEESSRESSERE REESE ESE FORESEES SSESEOSEEEFSE SEES FESS S 
C ACTION = FNUNOL CLESTING AND RENUMBERING SUPRESSED) 
CERERESTEREEREREESEERER SC EERESERE SEER ES FERSOEEERE SSE FESEORESSS ESET ESO FERS E 
600 WRITE( MASTER) FITEF 
640 J=LsLENTOT 
MRITE (MASTER) NEECOS)s JASFWKED) pSTEMOAL ID  ITESTT (SDs (PHL Jot} l= 
11519) 
649 CONTINUE 
c 
Ceeees EQEMAT™, 
C 


1) FORMATE THL »* SUD IATE MONE ASSUME I. ccc eNEW MATA FOLLOWSS®"*,///) 
Pas) FOOMATUABs 2X 221%52010) 
125 TOREMATISELO.3) 
$75 FARMAT ES Z"% CC CUPRECTED NATA Frou LINE ",13,°)9%) 
aro FTRMATIS UX L3¢ Xela UXei Fe 2 Xe Se LX ep LVF 9. S— 2XdD 
24 ENRMATEOS/y™ -CO LEUE S139" DELFTEN FROM DATA SET®) 
FUPMATO/e,"% CF 4 NEW LING WILL BF ADDED PREC FOING LINE ".14) 
165 FLRMAT IZ” CC NFW LENE %,93,% ADNFED TO DATA SET *) 
1% TARMATULHL » @ —€C MeAGeMD) TS AN ILLEGAL UPDATE FUNCTION®, /> 
= 80) 
42% TRRMAT(OZZ® ALL CHANGES PRECFOING APNVE FUNCTION WERE APPLIED*, 
1/slHL?) 
46) TORMATOLHL,* -€C ,935% 9) FS AN INVALID LINE NUMBER. NT CHANGES 
LAPPLIFD TN DATA. EXECUTION TERMINAT I ON Se ° eo “,/,1H1) 
5OL FORMATOLAL) 
BELG FARMAT(S LX, TB LX eV SolXel Se2Kol GelXelDlFIe3e2X)) 
420 FORMAT(IX,*PFC AGS AGE FEST #OONLUS STRAIN AT STRAIN AT * 
LAX UM STRAIN STRAING) 


[8 ® a a a  o nnn ernr nn 


a ET an 


540 TUPMATC TX, "NN. TEMP TEMP BREAK MAXUM" 

—STRESS FNFEOGY ENDURANT FM) 
550 FORMAT (45 Xe PSTRESSM eg LEX GS MNENSITY™) 

560 FORMAT (LX, §~---------- = +--+ + wenn 2 een een een 
ce ae ee ae ae ee a ne oe ee ") 

5h5 FORMAT (5xe"(4K) CFD CFD PSTD ae) ae) (Ps 
-1) (TAN) co) 

S70 EFNRMAT(® sseresseseessrsnt res eters See Se east 2 eee es sete stessessss 
leeer tas tetas Ste SSS te Se ees SS aeE Ss SSeS ST SESS eRe SSesr tse rss sets ssss 
Qreszess=) 


S70 FORMAT CLHE/LX,"ErPURIFR TRANSFORM INFRARED SPECTROSCOPY —- E41N PtH 
“YSICAL PROPERTIES MASTER TAPE GENERATOR® 
27UKe MCAFFPLZPCN — THINKOLSHUNTSVELLED®//1X,"PPNPERTIFS FOF SALIN P 
AROPFLLANT *,9X%,29AG) 
710 FORMAT (1HO," (COCL ACTION CARD NUMAER —*%913," 1S AN OUT OF SFOUERC 
TE yA, COMMAND. JOS TS TERMINATED. cco) 


Cc 
650) CONTINUS 
REWIND MASTF? 
PEWIND G 
PFTURNY 
FND 


LE TVET ERT aT, 
oe 


APPENDIX B 


PROGRAM FLOWCHARTS 


izslisee 
Card 10 


TABLE OF CONTENTS AND REFFRENCES 
PAGE se0x MANE 


FORTRAN NOOULE 


£490 ~ MAIN PROGRAR 


CHART VETLE — INTRODUCTORY COMMENTS 


CHART TITLE ~ PROCEOURES 


tooo00eo) 
teeo0n2) 
4000084) 
#000090) 
(0000035 
(0001023 
000104) 
(000108) 
1000106) 
4000112) 
(ooei19) 
(0001249 
1000128) 
FOOo1 sks 
1000133) 
(000135) 
(000136) 
1000137) 
t000142) 
(oo0150) 
1000153) 
1000182) 
(000184) 
toe@eo1n7) 
(ooeie6) 
(ooaz00) 
4000202) 
(000205) 
(000200) 
(ood. a: 
(800712) 
(0002201 
40002725) 
1060234) 
(000? +5) 
(002379 
1900299) 
(000244) 
(900248) 
1000232) 
(000254) 
cooozse) 
1000240) 
(e002se) 
(0022625 
(oo0267) 
(0002704) 
(ooe7ee) 
(e002 76) 
(900782) 
10902939 
toonz97) 
¢e00300) 
1000709) 
0000333) 
4000315) 
toeos1e) 
1000326) 
tows 201 
1000330) 
10003400 
(9803461 
(0003434 
(002356) 
$0036) 
18903579 
toon7en 
1o0e0361) 
(900349) 
(eaom7t) 


2006 
2.67 
2.00 
rrrty 
2.14 
2.23 
2026 
7.28 
2220 
2.0% 
2044 
74? 
2.48 
3.02 
3.03 
3006 
3008 
3210 
3oh8 
3.28 
4.08 
Hol? 
4.13 
4.17 
4023 
4e% 
4.26 
4.29 
oeo31 
3.¢1 
5.02 
5.03 
$07 
514 
5-15 
Sol? 
5.09 
Se2t 
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